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APPLICATION OF THE METHOD OP COORDINATE PERTURBATION 
TO UNSTEADY DUCT FLOW^ 

By Seymour C. Himmel 


SUMMARY 

The method of coordinate perturbation is applied to the unsteady flow 
of a compressible fluid in ducts of variable cross section. Solutions, 
in the form of perturbation series, are obtained for unsteady flows in 
ducts for which the logarithmic derivative of area variation with respect 
to the space coordinate is a function of the "smallness" parameter of the 
perturbation series. 

This technique is applied to the problem of the interaction of a 
disturbance and a shock wave in a diffuser flow. It is found that, for a 
special choice of the function describing the disturbance, the path of the 
shock wave can be expressed in closed form to first order. The method is 
then applied to the determination of the flow field behind a shock wave 
moving on a prescribed path in the x,t-plane. Perturbation- series solu- 
tions for quite general shock paths are developed. 

The perturbation- series solutions are compared with the more exact 
solutions obtained by the application of the method of characteristics. 

The approximate solutions are shown to be in reasonably accurate agreement 
with the solutions obtained by the method of characteristics. 


I. INTRODUCTION 

Problems involving the unsteady flow of a compressible fluid in ducts 
of variable cross section are frequently encountered in the study of 
nonsteady- state operation of air-breathing propulsion systems. At super- 
sonic flight speeds such problems are often complicated by the presence of 


*4rhe information presented herein constitutes the major part of a 
thesis that was offered in partial fulfillment of the requirements for the 
degree of doctor of philosophy. Case Institute of Technology, Cleveland, 
Ohio, June 1958. 
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shock waves in the flow. A typical example of a situation involving both 
an unsteady duct flow and a shock wave is that which arises when the 
throttle setting of a turbojet or a ramjet engine operating at supersonic 
flight speed is changed. As a result of the change of"- throttle setting, 
a disturbance is generated in the subsonic flow region of the diffuser 
that propagates upstream in the diffuser and interact^ with the shock 
wave. The interaction of the disturbance and the shock wave causes the 
latter to change its position in the diffuser and its strength. A knowl- 
edge of the history of the shock motion and- of the flow variables within 
the diffuser during such transient operation is of great value in the 
study of engine dynamics and the design of engine control systems. 

The difficulties encountered in attempting the solution of the par- 
tial differential equations describing the flow of a compressible fluid 
are well known. The usual procedure in solving problems governed by these 
equations is to~reduce them to a more manageable form by omitting terms 
whose effects are of small magnitude for the problem under investigation. 
In the case of ducts of variable cross section, such simplifications lead 
to the concept of quasi- one- dimensional flow; commonly employed in steady- 
flow theory. In the quasi-one- dimensional approximation, it is assumed 
that the cross-sectional area of the duct-caries slowly with distance 
measured along the axis of the duct. Under these conditions the velocity 
of the fluid is assumed to have the direction of the duct axis and all 
flow .variables are assumed to be uniform over any duct cross section. 

For unsteady quasi'-^ne-dimensional flows one' has to deal, therefore, with 
a single space coordinate and the time as the independent variables of 
the problem. Even with these simplifications the equations cannot, in 
general, be solved analytically. 

In the main, unsteady quasi- one-dimensional flow problems have been 
treated by the method of characteristics (cf . ref. 1 and the extensive 
bibliography therein) . Solutions are obtained by numerijcal methods or by 
a combination of graphical and numerical methods . In ^tther case much 
labor is involved and only the answer to an individual problem is ob- 
tained. A reasonably accurate approximate analytical method for treating 
unsteady quasl-one-dimensional flow problems is therefore desirable. 

Among those who have studied unsteady duct flow by analytical methods 
are ifentrowitz (ref. 2) and R. E. Meyer (ref. 3). Kantrowitz studied the 
foimiation and the stability of shock waves in duct flows by linearizing 
the equations of motion. He was able to demonstrate the_ instability of 
shock- free diffuser flows by this method. For diffuser .flows in which a 
shock wave is a part of the equilibrium flow, he was able to demonstrate 
the stability of the position of the shock.' Because he .was primarily 
concerned with stability considerations, Kantrowitz considered disturb- 
ances in the form of pulses, and his discussion of the interaction of a 
shock wave and a disturbance centers about conditions in the immediate 
vicinity of the shock. 
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Meyer treated waves of finite amplitude in ducts by considering the 
relations among the derivatives of the Riemann invariants along the 
characteristics. In this manner he developed a first-order theory for 
advancing and receding wave fronts in shock-free duct flows. By means of 
this theory^ he demonstrated the inherent instability of shock-free dif- 
fuser flows and was able to give a partial solution for the interaction 
of two wave fronts. 

In what follows, an approximate analytic method for treating the un- 
steady quasi-one-dimensional flow of a perfect fluid is developed. This 
method differs from those previously considered in that it is based on 
the method of coordinate perturbations (ref. 4). Solutions, in the form 
of perturbation series, are obtained for unsteady flows in ducts for which 
the logarithmic derivative of area variation with respect to the space 
coordinate is a function of the "smallness" parameter of the perturbation 
series . These perturbation series have as independent variables the 
characteristic parameters of the hyperbolic differential equations gov- 
erning the flow. This technique is applied to the problem of the inter- 
action of a disturbance and a shock wave in a diffuser flow. It is found 
that, for a special choice of the function describing the disturbance, 
the path of the shock wave can be expressed in closed form to first order. 
The method is then applied to the determination of the flow field behind 
a shock wave moving on a prescribed path in the x,t-plane. Perturbation- 
series solutions for quite general shock paths are developed. 

The solutions obtained by the approximate analytic method are com- 
pared with the more exact solutions obtained by the application of the 
method of characteristics using a finite-difference technique. The per- 
turbation series solutions are shown to be in reasonably accurate agree- 
ment with the solutions obtained by the method of characteristics. 

The author is very pleased to be able to take this opportunity of 
acknowledging his indebtedness to Professor G. Kuerti for his guidance, 
encouragement, and many valuable suggestions and criticisms throughout 
the preparation of this thesis. 


II. APPLICATION OP THE MESCHOD OP COORDINATE PERTURBATION 
TO QUASI- ONE-DIMENSIONAL PLOW 

In the usual perturbation theories for supersonic flow and wave 
propagation, the solutions are represented by perturbation series in 
which the flow variables are given as functions of the physical coordi- 
nates of the problem. Por example, for a small deviation from a uniform 
two-dimensional steady flow, the velocity components u and v are 
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given by expressions such as 

u = U + eu^^^(x,y) + e^u^^^(x,y) + • • •, 

( 2 . 1 ) 

V = ev^^^(x,y) + e^v^^^(x,y) + • • • 

(All symbols are defined in the appendix.) In these eq.ua tions, U is the 

undisturbed flow velocity and u^^^(x^y) is a deviation (in the x direc- 
tion) from the uniform flow. The magnitude of the deviation is governed 
by a "smallness” parameter e. For many problems this_method is adequate. 
In some cases, however, such solutions prove_to be inadequate; for exam- 
ple, in the Prandtl-Meyer expansion in two-dimensional steady flow (cf. 
section 2 of ref. 4). 

For hyperbolic differential equations the characteristic parameters 
are the natural Independent variables. In a perturbation theory based on 
the characteristic form of the differential equations, lit may be e:!^ected 
that the difficulties encountered in the more usual perturbation method 
may be avoided. Such a theory was developed by C. C. Lin (ref. 4) for 
quasi-linear systems in two independent variahles; it is based on ideas 
inplied in the work of E. E. Meyer (ref. 3) and suggested by K. 0. 
Friedrichs (refs. 5 and 6). It is referred to as the method of "coordi- 
nate perturbation.’’ 

In this method characteristic parameters a and are introduced, 
and the physical coordinates (say, x and y in two-dimensional steady 
supersonic flow) as weiJL as the flow variables u, v ate expressed in 
terms of the parameters. The perturbation solution thus appears in the 
fom 


u » u^^^(a,3) + eu^^^(a,^) + e^u^^^(a,3) + 
V = v^®^(a,p) ^-“ev^^)(a,^) + e^v(^)(a,p) + 

X = x^*^^(a,p) + €X^^^(a,p) + €^x^^^(a,p) + 

y = + ey^^^(a,p) + e^^^^(a,p) + 





( 2 . 2 ) 


• J 


Parametric, representations ^ such as that given above^ have disadvantages 
associated with the inapping onto the characteristic plane.. First, bounda- 
ry conditions are normally specified in the physical plane^ and it Is 
usually difficult to impose the boundary conditions in the plane of 
characteristic parameters. Second, a single-valued solution in the plane 
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of characteristics could become multiple- valued in the physical plane. 

However^ if^ as in the case considered^ the flov is a perturbation from a 

uniform state, the freedom in the choice of the characteristic parameters 
makes it possible to avoid multiple- valued mappings and to specify the 
boundary conditions easily. 

In the paper in vhich the perturbation theory vas developed, Lin 
treated problems of two-dimensional steady supersonic flow and emphasized 

the proper choice of the mapping, that is, of the parameters a an^ 

He was able to give a convergence proof for this case. In her thesis 

(ref. 7 ), P. Fox applied the method to the propagation of plane ("strictly” 
one-dimensional), cylindrical, and spherical waves. She was able to give 

a convergence proof for the case of plane waves. In what follows, the 

method of coordinate perturbation is applied to quasi-one-dimensionai 

flow. 


A. Equations of Motion 

The equations governing the homentropic, quasi- one-dimensional flow 
of a compressible fluid are (ref. 8): 



where A’ = A'(x') is the cross-sectional area of the duct at the station 
x’ . These equations are in terms of dimensional quantities. In what 
follows, it will be convenient to use dimensionless quantities for the 
variables given by: 





( 2 . 4 ) 


where the asterisk refers to the critical condition 'of the basic steady- 
flow problem; that is, where the local particle speed equals the local 
speed of sound. Transformation of equations (2.3) to equations involving 
the dimensionless quantities does not change the form of the equations 
but merely affects the boundary conditions. 
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When the dimensionless equations (2.3) are -written In canonical form, 
they become; 


= (u + c)t^, 
Xp = (u - c)tp^ 








- ac 


Ax • 

p = uo -j- tp,^ 


(2.5) 


vhere 


a 



( 2 . 6 ) 


and a and p are the characteristic parameters- ^us, (2-3) is re- 
placed by four differential equations for four unknowns in such a way that 
each equation involves only differentiation with respect to one character- 
istic parameter. 

While the pair (2.3) is irreducible in the sense of reference 6 ^d 
thus Riemann's method of integration can no longer be applied^ equations 
(2.3) or the equivalent system (2-5) can^ of course, -be solved by the' 
application of the method of characteristics, and many solutions for par- 
ticular problems have been calculated in this manner T(cf. for example, 
ref-, l) . It is advantageous, however, to develop an'_approximate analytical 
solution for perturbation problems governe_^ by these equations. Such a 
solution is obtained by the application of the method of coordinate 
per turba t i on . 
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B. Perturbation Solution 

A solution of equations (2.5) is sought in the form of pertu rbation 
series given by 


X = x^^^(a,p) + ex^^^(a,^) + e^x^^\a,3) + 
t = + et^^^(a,p) + e^t^^\a,p) + 
u = u^'^^Ca^p) + eu^^^(a,p) + e^u^^^(a,p) + 
c = c^^^(a,p) + ec^^^(a,p) + €^c^^^(a,p) + 




(2.7) 


The solution is, of course, a function of the duct area variation A(x) . 
In order to render the initial investigation as simple as possible, the 
duct area "was assumed to vary according to the relation 

A = (2.8) 


where b is a dimensionless constant and L is the dimensionless loca- 
tion of the duct throat (cf. (2.4)). For this choice of area distribution, 
equations (2.5) become ^ 


x^ = (u + c)t^^, 
Xp = (u - c)tp, 

+ ac = ucbt„, 
a (X ^ 


> 


(2.9) 
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Substitution of equations (2.7) into equations (2.9) yields 



( 2 . 10 ) 

Upon expansion of these equations and equating the coefficients of like 
powers of — e, there is obtained 
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Typical sets of these equations are: 
Zeroth order: 


(0) 

r (0) 


= (u 

xp 

= Cu^°) 

(0) 

+ c=Co) 

a 

a 

Jo) 

Up 

1 

Q 

O 

o 


. ^Co). (0) 

+ C ; 




a 


( 2 . 12 ) 


First order: 




_ = (u(°) . ^ Cut") * 

CCf CL 




X, 


( 1 ) = 

0> 

( 1 ) 


p 


. (0) (0) (1) (1) 

= (,u - c }tp + (u - c 


'a ■> 

(l)^.(o) 


)t 


P 


uC") . crc(l) = + ,(1),(0)),(0)1 I 

a a' ^ 

! 


> (2.13) 


- -h 

Up - acp = -b 


(0) (0) (1) , (0) (1) ^ (1) (o)u(o) 

uctp+^,uc +UC ;tp 


It is of interest to note that, whereas the zeroth-order equations 
are quasi-linear (they are identical in structure with (2.9))^ the higher 
order equations are linear and homogeneous in the unknowns of the same 
order; that is, the coefficients of the unknowns are functions of the 
solutions of the lower order equations. 

We wish to consider small unsteady perturbations of a steady flow in 
a duct; thus, the unperturbed flow should be a steady duct flow. This 
unperturbed flow is represented by the zeroth-order solution. The steady 
flow in a quasi -one- dimensional duct is cong^letely determined by the con- 
tinuity and Bernoulli equations in dimensionless form: 






u 


1 - 


u 


a + 1 


a/2 


2 2 

u + ac = a + 1. 


(2.14) 

(2.15) 
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These equations must be considered as time-^dependent integrals of the 
nondimensional form of equations (2.3) and^ since equations (2.5) axe 
equivalent to equations (2.3), they are solutions of these also. The 
zeroth-order equations (2.12) are the same.as equation's (2.5) with the 
area variation given by equation (2.8) used to evaluate the term i^/A. 

Thus, the solution of the zeroth-order equa,tlons is given by 


b(x^°^-L) 






a/2 




(2.16) 


^ ^^( 0)2 ^ 


a + 1. 


(2.17) 


In order to proceed to the solution of- the higher order equations, 
the zeroth-order solutions must be expressed as functions of a and p. 
To determine the form taken by these functf^s, an exaii^ple was chosen 
with 0 = 1 (t = 3) . This choice of y permits a simple algebraic form 

of the solution of equations (2.16) and (2.17) for u(®) and cC®) as 

functions of x^®^ . For subsonic flow downstream of the throat, these 
solutions are: 


JO) = . [x . 

(2.18) 

■ 

Joi-TlH- '\/l - 

L— 

(2.19) 


from which it can be shown that 

- 

— 

JO) . JO) , V3 [l - e^Cx'°)-L)]V3^ 

(2.20) 


JO) . e(0) = - Vs [l + 

(2.21) 


The zeroth-order characteristics are determined from the first 

two 


of equations (2.12). On a characteristic 6 = constant, we have 

.. .. 

— 


therefore : 


dx 


( 0 ) 


b(x^°^-L)i 


T/2 


= sjz dt 


( 0 ) 


( 2 . 22 ) 
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Setting 




(2.23) 


whence 


equation (2.22) becomes 


= ^b 


T- = -y/2 ■fa <it^°^. 

vl/2 ^ 


(2.24) 


5(1 - 5) 

Integrating this along the P = constant characteristic, we obtain 


(2.25) 


In 


- 1 


ViTT + 1 


= Vs bt^°^ + 


(2.26) 


At this stage the freedom of selection of the characteristic parameters 
mentioned previously enters. This means that the equations are invariant 
if one replaces a and p by a’ = f(a) and P' = g(P). It can be 
easily demonstrated that under these circumstances we may stipulate that 
a given curve of the flow plane (i.e., the x,t- plane) should be the image 
of a given curve of the a, P-plane. Now it is desired to have the 

zeroth-order solution represent the initial steady flow in the channelj 

that is, along the line t = 0, for example. Therefore, we select the 
characteristic parameters such that, on t^O, x=a=P (i.e., the x-axis 
is the image of a - p = O) . For this choice of the parameters, equation 
(2.26) becomes 


In 


Vl^ - 1 


VT^ + 1 


- In 


Vl - C(P) - 1 


VI - 5(.3} + 1 


= -y/2 bt 


( 0 ) 


(2.27) 


Similarly, on an a = constant characteristic we have 


In 





Vl +5(0,) - 1 

- In 

VTTl"- 1 

_yi + + 1 

-^1 + ^ -t- 1 


= -y/2 


( 0 ) 


(2.28) 
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Eliminating from the last two equationa and solving for $ 

obtain 


where 


and 




4 1 + I 

CtI - 1)^ 




s = 


Vl - ■ 1 


„-a£IS^La 

1 

From these equations ve obtain for the physical 


coordinates : 


X 


( 0 ) _ 



-V^ , 

(ti - if 


Vs bt^®^ = In 



- l)^ - .^(n + |)V^ - (ti - _il . 

'^(n - 1)2 - 4 (ti + |)V^ + (ti - |)_ 


By substituting equation (2.29) into equations (2.18) and (2.19), 
particle velocity and speed of sound are given by 

u<°> = - nrHr[('' - 5)^ - V*! - s)^ 

yCn - s)* -neCti + 6)^ 


we 


(2.29) 


(2.30) 


(2.31) 


(2.32) 


(2.33) 

the 


(2.34) 


(2.35) 


6£?9 
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These rather formidable equations give the parameters of a steady- 
duct flo-w (l.e., the zeroth-order solution) in terms of the characteris- 
tic parameters a and 3. The higher order approximations are given by 
linear equations where the coefficients are functions of the lower order 
solutions. For example, the first-order equations (2.13) have coeffi- 
cients derivable from equations (2.32) to (2.35) . The higher order equa- 
tions are, howe-\rer, complicated and must be solved simultaneously. No 
simple analytical technique seems applicable - perhaps a numerical tech- 
nique would be feasible. 


To illustrate the nature of the zeroth-order solution for the area 
variation selected, an example -was computed. The flow situation is J.1- 
lustrated in the sketch. 



We ha-ve a duct with area increasing to the left . . The throat of the duc t 

is located at x^®^ = L, and the flow is from right to left. The coeffi- 
cient b -vrais chosen as 1, and L -was taken as 2. Typical characteris- 
tics of the zeroth-order solution for subsonic flow do-wnstream of the 
throat are shown in figure (2.1). The characteristics ^ = constant are 
essentially straight at the left, that is, the region of large areas, and 
have increasing curvature as they approach the throat, where they have a 
vertical tangent. The characteristics a = constant are es sen tially 
straight throughout the duct. The variations of particle and sonic speed 
and Mach number along the duct are shown in figure (2.2) . 


At this point it is of value to digress from the par ticu lar problem 
at hand and to note an important feature of the method of coordinate 
perturbation that is implicit in the preceding development. In the case 
under consideration, we find that cc, and p niay, at the same time, be 
considered as the characteristic parameters of a zeroth-order s olution 
and of an "exact" solution (in the example, the steady duct flowT^ That^ 
is, the perturbation-series solution begins -with the exact characteris- 
tics of the underlying flow. The lines a = constant and p = constant 


are mapped by x^®) (a-,P) and t^°^ (a-, P) onto the characteristics in the 

x(0),t(0) -plane and by x(a, p) and t(a, p) onto the characteristics in the 

x,t- plane. Had it been possible to determine the solutions of the 

equations of higher order, these higher-order terms would not only have 
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changed the distribution of the flow variables along the zeroth-order 
characteristics in the physical plane but would also have changed the 
characteristics themselves. That is, the lines a = constant and 
P = constant would no longer be mapped by x(a, p) and t(a,p) onto the 
characteristics of the underlying flow in the x,t-plane but instead - 
would be mapped onto the approximate characteristics _pf the perturbed 
flow. In other words, the significance and the properties of the charac- 
teristics are retained in the method of coordinate perturbation. This is 
in marked contrast to the more usual perturbation theories in the physi- 
cal plane, in which the characteristics of the perturbed flow do not ap- 
pear explicitly. 


C. The "Epsilon Duct" m ■ 

A way of avoiding the difficulties involved in obtaining the higher 
order solutions for the preceding example presents itself when the duct 
area variation itself is considered as being of the same order as the 
"smallness" parameter e. In particular, the case A^/A = - € was 

studied^ such a duct shall be referred to as an "epsilon duct." Such a 
choice amounts to requiring that the duct area vary slowly, but places no 
restriction on the over-all area variation (i.e., over a long stretch of 
the duct). For this choice of area variation, equati^s (2.9) become 


x^ = (u + c)t^, > 

Xp = (u - c)tp. 

Up - acp = - €uctp.^ 


(2.36) 





O ctO «J 


NACA OM 1439 


15 


Substituting equations (2.7) into equations (2.36) and equating the coef- 
ficients of like powers of e yields 


OCi 


g . c^^^)t^^'^^ k=0,l,2 

= Oj 4°^ - = 0; 


(k) (k) 


k -1 f J 

= E E 

j=0 A=0 


(j- 7 ^) *Ck-l-j) 


J=0 A=0 _ 


; ]s?=l,2,3 


Typical sets of the approximate equations are: 
Zeroth order: 


(2.37) 


= (uCo) . 


xf ) = (u'°) - 

4 °^ " 


(2.38) 
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First order; 


Cu (X Cu 

41) = („C0) _ ,(0)„^1) ^ („(1) . ,(1)),(0)^ 


4^’ 


+ ac 


( 1 ) .(o)..(o).(o) 




a 


= c u 


'"a ^ 


ac 


(1) = . 


(2.39) 




p ^ ..p 

Second order: 

.(2) . (JO) ,.J 0 ))J 2 ) ^ (^(1) ^ ^(1)),(1) ^ ^J 2 ) ^ eC2))t<°), 




X. 


a 


a 


a 


= (u(°) - 




p 


(a) , , 42 ) . eC 0 )J°)tCl) . 


>(2.40) 


( 2 ) 

^P 


. 0 ^ 2 ) _ , 

- acp - - c 




J 


Third order: 


,( 3 ) . (Jo) ^ ,,(0))^.(3) ^ (Jl) ^ + (u(2) + o(2))tCi) + (u(=) + c(3))t(°), 

,(3) ^ (JO) . J0)j,(3) ^ (Jl) . Jl)j,(2) ^ (J2) . J2)j,(l) , (J3) . J3)j,|0)^ 

H- - ctO)utO)t(2) . (e(O)u(^) . . =(">u(l) H- c<2)JO))JO)^ 


> (2.41) 






(o)Jl) ^ Jl)JO)),(l) . (J0)J2)^ JDJl) ^ J2)J0)),(0). 


- (o'* 


The interesting thing about these equations is that the right hand 
sides of the equations resulting from the compatabillty conditions 
(i.e., the third and fourth equations of each set) involve functions of 
the lover order solutions onlyj that is^ if the latter are knovn^ of the 
independent variables. Thus^ in principle j. these equations can alvays be 
solved by quadratures. The equations resulting from the characteristic 
conditions (i.e., the first tvo equations of each set) are also alvays 
linear^ and the coefficient of the derivative of the'k-th order approxi- 
mation of t is alvays the sum or difference of the! zeroth- order solution 
for u and c. -■ I" 
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The form of the equations suggests that it vould he particularly 
simple to obtain solutions for problems where the Initial conditions are 
specified on a line x = constant j for in such cases the zeroth-order 
solutions for u and c would be constants. This renders the solutions 
for the physical coordinates relatively simple for an orders. 

To demonstrate such a solution, and in order to be able to check the 
values obtained against a known exact solution, it was decided to solve 
the problem wherein it is specified that on the section x = 0 the values 
of the particle and sonic speed are given by the constants tJ and C, 
respectively (i.e., a time- independent flow). The constants U and C 
are subject to the condition 

= or + 1, (2.42) 

that is, the Bernoulli equation. The solution generated under these con- 
ditions should be the steady-state duct flow. 

Equations (2.38) to (2.41) are to be solved, therefore, under the 
boundary conditions given above. The boundary conditions are to be 
satisfied in the following manner: Contrary to the choice in the previous 
section, now the line x = 0 is to be mapped on a = p such that 
a = p = t. Thus, we choose 


x^“^(a,a) = x^^^(p,p) = 0 for 

n ^ Oj "N 


t^^Ha^a) = a, t^°^(p,p) = p. 


(2.43) 

t^^^(a,a) = t^^^(p,p) = 0 for 

n> 0} J 



and 


u^°^(a,a) = u^®)(p,p) = U, 
c^®Ua.,a) = c^®^(p,p) = C, 

u^^\a,a) = u^^^(p,p) = c^^^(a,a) = c^^^(p,p) = 0 for 

The second pair of the zeroth-order set (eqs. (2.38)) gives 
particle and sonic velocities slnply 

u(®) = Uj cC®) = C. 


N 

> (2.44) 


n> Oj 
for the 
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Inserting these in the equations for the characteristics 
yields 


- (U + C)t(O) = 

x(0) - (u - C)t^°) = g?(a). 

Applying the boundary conditions, we obtain 


^(8) = - (U + C)p, 
g^'Ca) = - (U - C)a. 


Therefore, 


.CO) = 


.( 0 ) = s - a 2 


2C 


20 


The solutions of the higher order sets are obtained in a 
The solutions up to the third order are ' 


u(°) . a, 

- (a - ^) f, 

= (a - p)2 ^ (oo2 . u2), 

= (a - p)® ^ (oO^ - SU^)i 




> 




C^Q) = C, -N 

c(2) = . (e . > 

o(5) - - (a, - P)’ -2^ (20C^ - U®); 

24cf^ 


and integrating 


similar manner. 


(2.45a) 


(2.45b) 
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jr 


o> 

fcO 

cvj 

LO 


M 

O*- 

M 

to 

I 


.( 1 ) 


.( 0 ) 


= - (a - 3)' 


(a - p) 


2C ^ 




jjI 

8aC" 


C^( 2 a + 1 ) + U^, 


= - (a. - p)= 


_uL 

24oC 


(gC^ + 

gC" 




+ (g + i)c 


.2 

t 


> (2.45c) 


.( 3 ) 


= - (a - 3)^ I |j:K - N) g + (K + j 




vhere ^ 


+Co) _ o tj + c u - c 

^ “^20 " 20 ' 


tW = . (a - s)2 _2_ (oO^ + D^), 

8 gC 2 




t( 2 ) = . (a - 3 )" "" 


48 g^C^ 


(gC^ + U 2 )(gc 2 + 2 U^), 


> (2.45d) 


t^ 3 ) = _ (a - 3 )^ ^ (K - N) 


K - N = — ^ (gc2 + U2)(cc2 + 6U2), 
48 g^C^ 


J 


(2.46) 


2 2 r* — 

K + W = - (gC^+U^)(2g^c2+u2) + JL_ a(g+4)c2 _ (Sg+2)u2 .(2.47) 
ISg'^C^ 48 g^ L J 


As anticipated^ the solutions for the flow variables are independent 
of the time as is seen from the fact that c^ and x are functions of 
(a - P) alone^ and the possibility of the elimination of (a - p) from 
these equations makes u and c functions of .x alone. 

Although attempts were made, it was not possible to prove analyt i cally 
the convergence of the perturbation series. In any case^ however, it 
would have been necessary to resort to numerical examples in order to in- 
vestigate the accuracy of the representation in the various regions of x. 
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As an example^ the duct e = 0.1, L = 6 was chosen; that is, the duet 
area variation is given by A = exp^-0.l(x - 6)J . A plot of the duct area 

variation is given in figure (2.3). The cruss-sectional area varies in 
the ratio of approximately 1.81:1 from initial line (x = O) to throat. 

This duct was first investigated for the case of subsonic flow down- 
stream of the throat with o = 1 (r = 3) . With the flow coming from the 
right, we have, from the steady-flow equations (2.14) and (2.15), on 
X = 0 


U = -0.40504; C = 1.35497. 

Inserting these values in equations (2.45), we obtain the series solution 


-u = 0.40504 + 0.02744(a - 8) + 0.00085(a - ^)^ + 0.0000l(a - 
c = 1.35497 - 0.00820(a - p) - 0.0005S(a - p)^ - 0.00002(a - 
X = 0.61695(a - p) - 0.00633(a - p)^ - 0.00029(a - p)® - 
0.0000l(a - p)^, 

t-= 0.35053P + 0.64946a + 0.00552(a - p)^ + 0.00015(a - p)^ + 
0.000003(a - P)^. 



>(2.48) 


J 


The second- and third-order solutions are shown in' figures (2.4) and 
(2.5) by the data symbols. The exact steady duct flow is given by the 
solid curves. As can be seen in figure (2.4), the flow velocity is very 
well represented by the second-order solution up to about x =» 4. The 
third-order solution extends the accurate representation up to x = 5.5. 
Only in the Immediate vicinity of the throat (x = 6) does the departure 
from the exact solution become relatively large. The sonic speed is also 
well represented by the approximate solutions, as seen in figure (2.5). 
Again, large departures from the exact solution occur only near the throat. 


The variation of the space coordinate x with the" quantity (a - p) 
is shown in figure (2.6). Proceeding from the second-. to the third-order 
solution decreases the variation of x with (a - p) for the larger values 
of this quantity. That is, the higher degree terms of __(a - p) cause the 
X, (a - p) -relation to become more horizontal for values’ of the argument , 
above 7. For values of (a - p) less than 7, the second- and third-order 
solutions give the same value for x. It may be recalled that values of 
(a - p) which yield values of x greater than 6 are not pertinent to the 
solution . 
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A typical p = constant curve in the x^t-plane of the third-order 
solution is shoTO in figure (2.7) along with the exact characteristic. The 
case shown in that of p = 0^ and the points indicated by the data symbols 
are obtained by setting P = 0 in the last two of equations (2.48) and 
eliminating a. The exact characteristic is obtained by setting b = e 
in equation (2.27) and using the appropriate value of L. It is evident 
that the approximation to the actual characteristic is very good up to the 
vicinity of the throat, in which region the departures become Isirge. 


D. Factors Affecting Convergence and Accuracy 

From the preceding example, it appears that the approximate solution 
does converge to the duct flow solution. There are a number of factors 
that may affect the nature and rapidity of the convergence. Among these 
are the magnitude of €, the over-all area ratio of the duct, the loca- 
tion of the Initial line on which U and C are specified, and whether 
the flow is subsonic or supersonic. 

1. liyf ect of € . - The first step in this investigation was to deter- 
mine the effect of the magnitude of e on the solutions. For this pur- 
pose a set of approximate duct flow solutions were computed for different 
values of e. The values of duct length were so chosen that the product 
€L was constant. For the type of duct considered, that is, 

A = exp j^-e(x - L)J , this results in a set of ducts of the same over-all 

area ratio between x = 0 and x = L but with different rates of change 
of area. These examples were computed for the case of subsonic flpw down- 
stream of the throat with cr = 5 (r = 7/S) • The values of the parameters 
€ and L used are; 



For these cases we have, from the steady- flow equations, the boundary 
conditions on x = 0; 


U = -0.36841; C = 1.08298. 
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For these values the approximate solution is given by 

-u = 0.36841 + €(0.19949) (a - P) + €^(0.05276)(a - ^)^ + 
€^(0.00862) (a - 

c = 1.08298 - €(0.0l357)(a - 8) - €^(0 .00735) (a - 8)^ - 
e^(0.00262)(a - 8)^^ 

X = 0.47883(a - 8) - e(0.03772)(a - 8)^ - e^(0.01376)(a 
€^(0.00038) (a - 8)^. 




X2.49) 



The results of the computations for these examples are shown in fig- 
ures (2.8) to (2.10) as plots of u and c as functions of x. For «n 
the -values of e, the approximate solutions converge on the exact solu- 
tion. The accuracy of the representation is quite good^ the third-order 
solution being quite accurate for about 80 percent of the duct length. 
Moreo-ver, comparison of the three figures sho-ws that the accuracy of the 
representation at corresponding points (i.e., equal area ratios) in the 
ducts is the same for all values of e. This leads to the conclusion 
that the accuracy of the approximate solution is not a function of e 
alone. --- ^ 


2. Effect of over-all area ratio . - The over-all area ratio of an 
"epsilon duct” is a function of the product eL. To “investigate the ef- 
fect of the over-all area ratio on the solutions, two additional examples 
■were computed for an € of 0.2. For subsonic flow the values of L 
selected and the corresponding Initial values of the. flow, -variables on 
X = 0 are; 


L 

u 

C 

5.0 

7.5 

-0.23889 

-0.14270 

1.09022 

1.09365 
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The solution for L = 5 is given by 

-u = 0.23889 + e(0.13022)(a - p) + e^(0.03515) (a - p)^ + 

€^(0.00614) (a - p)^, 

c = 1.09022 - e(0.0057l)(a - p) - €^(0.003H)(a - p)^ - 
€^(0.00113) (a - p)^, 

X = 0.51894(a - p) - €(0.01576)(a - p)^ - £^(0.00575)(a 
(0.00157) (a - p)^. 

-For L = 7.5 the solution is 

-u = 0.14270 + e(0.07803)(a - p) + 6^(0.02126) (a - p)^ + 

6^(0.00382) (a - p)^, 

c = 1.09365 - 6(0.00204)(a - p) - 6^(0.001H) (a - p)2 - 
6^(0.00041) (a - p)2, 

X = 0.53752(a - p) - 6(0.0056l)(a - p)^ - 6^(0.00205)(a 
6^(0.00056) (a - P)^. 

The results of these computations are shown in figures (2.11) and 
(2.12), where the particle and sonic velocities are plotted as functions 
of X. Again, the approximate solutions converge on the exact solution. 
Comparison of these solutions as well as that shown in figure (2.9), which 
is for L = 3 and e = 0.2, indicates that the greater the length of the 
duct (i.e., the over-all area ratio from x = 0 to x = L), the smaller 
is the value of x/L at which a given departure from the exact solution 
occurs. For example, the coordinate x at which the third- order solution 
for u has an error of 0.01 is given in the following table for the three 
ducts considered: 




>(2.51) 


- P)^ - 






>(2.50) - 


^)^ - 


L 

3c 

x/L 


(tij) exact 

3.0 

2.64 

0.88 

1.70 

0.760 

5.0 

3.80 

0.76 

2.14 

0.573 

7.5 

4.90 

0.65 

2.66 

0.405 




24 


mCA TM 1439 


Although a given error in u occurs at smaller values of x/L as the 
duct length increases, it should be noted that the ratio of the area at • 

X = 0 to that at x Increases as the duct length increases. A similar 
statement holds when the comparison is made on the basis of a given per- 
centage error in u. 

For a given value of €, the value of L determines not only the 
over-all areh ratio of the duct but also the values of the flow variables 
at X = 0, since the steady flow is sonic at x = L. The larger the value 
of L, the smaller is the value of U on x = 0. A concomitant of small 
values of U is a decreased rate of change of the flow variables with x 
(or area ratio) in the vicinity of x = 0 (cf. figs. (2.9) and (2.12)). 

It is this latter property of subsonic flow that results in the accurate 
representation of the flow over larger area ratios for the "epsilon ducts" 
with larger L. 

3. Effect of the location of the initial line . - To this point, only 
the case wherein the initial conditions have been specified at the largest 
area of the duct, that is, at x = 0, has been considered. If the initial 
conditions are specified at the throat, that is, at ~x = L, where the 
rates of change of the flow variables with distance are greatest, the 
effect of this factor can be observed. For this purpose, equations (2.36) 
were solved with the conditions specified at x = L. The solutions for 
this case are the same as those given by equations (2.45) with the ex- 
ception of the zeroth- order term for x. This is replaced by 

x(0) = L - (a - p) (2.52) 


For this case we have at the boundary x = L the initial values of 
U = -1^ C = 1. The steady- flow solution for cr = 5 "is represented by 
the following equations: 

-u = 1 + €(0.5)(a - P) + e^(0.l)(a - p)^, 
c = 1 - e(0.l)(a - P) - e^(0.05)(a - p)^ - e^(0.015)(a - p)^, 

X = L - €(0.3)(a - p)2 - e2(0.1l)(a - p)^ - 

€^(0.03075)(a - p)^. J 

These equations represent the supersonic branch as well as the sub- 
sonic branch. To obtain the subsonic branch, the argument of the series 

(a - p) must be taken as negative to make < 1. An example of the sub 

sonic case was computed with e = 0.1 and L = 6, and the- results are 
shown in figure (2.13). The representations of the particle and sonic 


M2.53) 


Rrv7r* 
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velocities given by the approximate solutions are seen to alternate above 
and belov the exact values. The approximate solution does, hoover, ap- 
/ pear to be converging to the exact solution as higher order terms are 

used. In the neighborhood ot the throat the approximate solutions are 
good representations of the flov, despite the rapid variation of the flov 
parameters vith x. Although the solutions do appear to converge to the 
exact solution, the convergence is not as rapid as that obtained for the 
corresponding case with conditions specified at x = 0 (fig. 2.8). This 
is especially noijiceable for the sonic speed (fig. 2.13b), where the 
second-order solution reaches a maximum at about x = 4 and then de- 
creases with increasing x. 

Some insist into the reasons for this behavior can be obtained from 
figure (2.14), which shows the variation of x as a function of (a - p) 
for this example. The second-order representation of x has a minimum 
in the vicinity of (a - p) = -17, whereas the third-order representation 
exhibits no such behavior up to this value of the argument. This change 
in the nature of the relationship between x and (a - p) is a result of 
^ the fact that, for negative values of the argument, the expression for x 

is in effect an alternating series. The coefficients of the series are 
of sufficient magnitude to have large effects on the representation of x 
for the larger values of the argument. It should be noted that for nega- 
- tive values of the argument the series for u and c also have alter- 

nating signs. 

Thus, for subsonic flow, it can be concluded that specifying the 
initial conditions at x = L, where the rates of change of the flow varia- 
bles with X are high, results in less rapid convergence of the solutions 
than for the corresponding case where -the initial conditions are given on 

X = 0. * 


4. Supersonic flov . - Thus far, only subsonic flows have been con- 
sidered. We now examine the character of the perturbation solutions for 
supersonic flow. First, -consider the case with initial conditions given 
at the throat. For positive values of (a - S), equations (2.53) yield the 
supersonic branch of the duct flow. In this case the series do not have 
alternating signs. The results for € = 0.1 and L = 6 are shown in 
figure (2.15). In this case the approximations converge quite rapidly to 
the exact solution, and the accuracy of the third-order solution is quite 
good for the entire length of the duct. Whereas in the subsonic case 
specified on x = L convergence of the representation of c was poorer 
than that of u, in the supersonic case the opposite is true. In fact, the 
third-order representation of c is excellent for the entire length con- 
sidered. One possible reason for the fact that the representation of u 
does not appear to be as accurate as that for c is that the coefficient 
of the term, of third degree in (a - p) in the equation for u is zero ' 
in this case. 
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There remains one further case to examine « This is supersonic flov 
with conditions specified at x «= 0. For this case we have, for the same 
duct as in the last two examples, on x = 0, U = -1.67289 and 
C = 0.00019. For these values the flow is represented by 


-u = 1.67289 + e(0.6693l)(a - ^) + e^(0.01685)(a - 3)2 - 
6^(0. 06018) (a - 3)^, 

c = 0.80019 - e(0.27986)(a - 3) - 62(0.H197)(a - 3)^ - 
6^(0.01681) (a - 3)^, 

X = -1.34859(a- 3) - e(l.07539)(a - 3)^ - e^(0.43969) (a - 3)® 
6^(0.16413) (a - 3)^. 


(2.54) 


For values of x greater than zero, (a - 3 ) fflust be taken as nega- 
tive. This again leads to a series with alternating signs for the rep- 
resentation of X. The results of the evaluation of equations (2.54) are 
shown in figure (2.16). For this case, as for subsonic flow specified at 
X = L, the approximate solutions of increasing order alternate above and 
below the exact solution. Again, the third-order solution is a good rep- 
resentation of the flow variables over most of the duct- length. Depar- 
tures of appreciable magnitude for the third-order solution occur in the 
vicinity of x * 5.3. It is to be noted that both the first- and third- 
order solutions for the flow variables fold back on themselves, the former 
at about x = 4.2 and the latter at about x = 5.5. This is again the 
result of the nature of the x, (a - 3)~relation shown in figure (2.17). 
The first-order representation reaches a maximum of 4.^ at (a - 3 ) ® -6.25, 
and the third-order representation has a maximum of 5.55 at (a - 3 ) = 

-9.0. These maxima result from the strong influence of: the term of 
highest degree of (a - 3 ) in each case. These are the even-powered terms 
that, for negative values of the argument, serve to reduce the rate of 
change of x with (a - 3) and thus produce the maxima. , 




5. Recapitulation . - From the preceding examination of the factors 
affecting the convergence and accuracy of the perturbation-series solution 
for steady duct flow, it was established that for both supersonic and 
subsonic flows the series converge to the exact solution. The extent of 
accurate representation of the flow depends on the over-all area ratio of 
the channel and thus on the values of the flow variable's on the initial 

line. For subsonic flow the extent of accurate representation was greatest 

when the Initial values were given at the large end of the duct, that is, 
at X = 0. For supersonic flow the opposite was true; that is, the extent ■ ■ ■ 
of accurate representation was greatest when the Initial conditions were 
given at the throat. However, the differences caused by the change of the 
data carrier in the supersonic case are not as great as for subsonic flow. 
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In both supersonic and subsonic flow the poorer representations occur 
when negative values of the argument of the series must be employed. 

It should be noted that the particular area function chosen, that 
is, the "epsilon duct," is not the only one for which such solutions are 
possible. The particular case used should be viewed as but one of a fam- 
ily of area functions that should be amenable to the method used. This 

family can be represented by A = where P(x) is a polynomial in 

X. For such area functions similar results should be obtainable. 


E. Unsteady Plow 

To this point, the perturb at ion- series solution h^ been investigated 
only for the case of constant conditions on the initial line. Such ini- 
tial conditions were chosen so that the solutions obtained could be easily 
compared with a known exact solution, the steady duct flow. We now turn 
to the case of unsteady quasi-one-dimensional flow. For perturbation- 
series solutions for unsteady flows in a duct, the specification of 
boundary conditions becomes complex, and more detailed consideration must 
be given to the nature of the initial data than was the case for steady 
flow. 


1. Boundary conditions . - In the usual perturbation problems, such 
as those for steady two-dimensional flow and one- dimensional unsteady 
flow, there is an underlying uniform state or flow that is constant both 
in space and time. For example, in two-dimensional steady flow the under- 
lying uniform flow is frequently chosen as a constant value of u, as in- 
dicated in equation (2.1). For one-dimensional wave propagation problems, 
the underlying flow is either a constant steady flow or a state of rest. 

In such cases the physical interpretation of a given set of boundary data 
is quite straightforward and no great difficulties arise. 

In contrast to these cases, the underlying flow in a quasi-one- 
dimenslonal problem is uniform, in the true sense of the word, only for 
one particular circumstance, that is, a state of rest. For all other 
cases the underlying "uniform" state is a steady duct flow in which, of 
course, the flow variables are not constant in space. For such problems 
an Implicit boundary condition is then: Along some line t = constant 

in the x,t-plane, a steady duct flow must exist. Alternatively, one could 
specify that, along a certain characteristic which separates the regions 
of steady and unsteady flow (corresponding to the "rest" characteristic 
of a disturbance advancing into a gas at rest) , the steady duct flow rela- 
tion holds . This requirement Introduces some difficulties in the specifi- 
cation of Initial data for the solution of problems by the perturbation 
method under consideration . . : 
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Consider tvo intersecting noncliaracteristic arcs J and K in the 
x^t^plane enclosing a region R (see sketch). 



As developed in section 28 of reference 6, if a direction is assigned 
to each of the two families of characteristics (in the sketch the direc- 
tion of increasing time is chosen)^ it can happen that the two character- 
istics issuing from- any point on K enter R^ whereas "only one charac- 
teristic from any pointnon J enters R. An arc such as K is called 
space- like^ and the arc J is called time-like. In the reference cited 
above^ it is established that^ for problems In which the boundary condi- 
tions are specified on both time- like and space-like arcs^ two data on 
the space-like arc and one datum on the time-like- arc determine an unique 
solution in the region between the arcs. In such problems^ a relation*, 
exists between the two dependent variables on the time- like arc. 

The form taken by the differential equations for the perturbation 
problem (eqs. (2.37)) is such that-it is most convenient to specify the 
initial data on a time-like arc. This was the case for the steady- flow 
examples presented earlier- in this section. In that case it was reasona- 
ble to specify both dependent variables on the time- like arc^ x = 0 or 
X = L, because the relation between the dependent variables was known and^ 
further, they were constant. In the more general case. of an unsteady 
disturbance of a steady flow, one cannot arbitrarily specify two data on 
a time-like arc as noted above. If such a specification were made, the 
region S, bounded by the characteristics issuing from point 0, would 
contain an unsteady flow. This would iTsply that a disturbance- originating 
at^oint 0 would influence the flow field to the right of 0 at an 
earlier time, a physically untenable condition. On the other hand, if it 
were attempted to specify a steady flow on the bounding C+ characteris- 
tic (the "rest” characteristic) issuing from_point 0 and a single datum 
on the time- like arc, the problem would become almost hopelessly complex 
because the initial data along the "rest" characteristic would have to be 
given by equations such as (2.30) to (2.33). 

There are, however, two ways of avoiding the difficulties noted above. 
The first is to consider the problem wherein two data are arbitrarily 
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specified on a time- like arc from a different viewpoint. The second is to 
seek the relation between the dependent variables on the time-like arc 
that will not produce an unsteady flow below the ’’rest" characteristic. 

•i 

2. The "inverse" problem . - ¥e consider first the case of an arbitrary 
disturbance of one or both of- the dependent variables specified on a time- 
like arc. Such a problem can be viewed in the following manner: We give 

■ ' on a time- like arc, say on the line x = Xq, a certain distribution of u 
and c, one or both functions of time. We then ask what occurred at 
another station at a previous time to cause the flow at x = xq to vary 
in the manner prescribed. For the flow to have been a steady duct flow 
prior to some time tQ on the boundary, it is necessary merely to specify 
that for t < to the values of u and c on the boundary are to be 

constant and satisfy the Bernoulli equation. This situation is shown in 
the sketch. 



Under these conditions there will be a steady duct flow below the C- 
characteristic through (xq^ to) . A physical example of such a problem 
would be that of inquiring how to vary the entrance conditions into a duct 
so that a given flow variation at the exit might result. 


Such problems are referred to as "inverse" problems in the sequel. 
It was found that the solution of such problems could be written in a 
quite general fashion. An example of such a solution follows. 


Consider the case of an "epsilon duct" originally containing a steady 
flow. Let a disturbance in sonic speed of order of magnitude € start 
at a specified time at a given cross section of the duct. For convenience 
let the cross section be that at x = 0 and let the disturbance begin at 
t = 0 at this cross section. The boundary conditions under which equa- 
tions (2.37) are to be solved are, therefore, on^ x = 0, 



(2.55) 
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vhere 


f (t) = 0, for t < Oj 
+ aC^ = a + 1. 


(2.56) 


These boundary conditions are to be satisfied in the following man- 
ner; The line x « 0 is to be mapped on a = p such that a ** p = t^ 
that is, according to eq.uations (2.43) and 

J°^(a,a) = u(°)(p,p) = U, ^ 

u^^^(a,a) = u^^^(p,p) =0, n > Oj 

c^*^Hct.,a) = o^°^(p,p) = C, 

o^^)(a,a) = f(a), c^^^(p,p) = f(p), 

c(“)(a,a) = c(“)(p,p) = 0, 




n > 1. 


(2.57) 




(Note the peculiar way in which the boundary conditions for this problem 
can be accommodated in the e- scheme.) 

The solutions, under these conditions, up to second order are: 

u(0)»U, 


- p) + U^^^f(p) - u|^^f(a), 

y(2) 3. u^2)(^ _ pj2 ^ u|^)(a, _ p)f(a) - - p)f(p) 

f(T)dTJ 


> (2.58a) 


J 


.(0^ - c 

C = C, 




- p) + f(p) + C^^^f(a), 

= - Cp^(a - p)^ + 0 ^^ (a - p)f(p) + 4^^(a - p)f(a) + 


>(2.58b) 
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= - 4°)(a - P), 

- P)^ + xl^^co - p)f(p) + X^^Ha - p)f(a) + 

4^^ f (T)<iT, 

x^2^ = - X^^^(a - p)^ + X^^\a - p)%(p) + X^2^(a - p)^f(a) + 
(t - p)f(T)dT + (a - T)f(T)dt + 

Xg^^f(p) f(T)dT + X^^^f(a) f(T)dT + 

X^^^ y^^jfCr^ ^ dT + X^^^Ca - p)[f(p^ ^ + 

X^Q^Ca - P) jf(a^ ^ + X^\a - p)f(a)fCp); 

^(O) ^ i(0)p _ 

- p)^ + - P)f(p) + 0?^^ka. - P)f(a) + 

f('r)dT, 

t(2) = . _ p)3 ^ gj(2)(^ _ p)2f(p) + tC 2)(^ _ p)2f(a) + 

^4^^ ■ P)^CT)dT + 

Tg^^fCp) f(T)dT + T^^^f(a) f(T)dT 

Tg^^Ca - P)jf(p^^ + 

T^2)(a,_. p)Jf(a^2 + t(2)(cc - p)f(a)f(p). 


>(2.58c) 


It + 


M2.58d) 
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In obtaining these solutions, the following identity has been used: 

S • •• T^rn'tP^ (y-£)^°~^^f(i)<is. (2.59) 


n 


n 


The coefficients appearing in equations (2.58) are functions of U, 


and 

are listed below; 





u(") 

UC 

2 

uC2) , 

. i (,0= - « 

2) 


0 


_ (u + 

c)(u - 

oC) 

2 

2 

^2 


80 


0 

U^2) . 


cKp-<- 

aO) 

3 

2 

3 


80 “ 




u<2) . 

_ P + - 





4 

4 


J 

^1 

U2 

2a 

0i . 

. ^ 
4a 



c(i) 

^2 

1 

=3 — 
- 2 

.(2) 
^2 ■ 

= |-(U 

- 0)^ 


c(i) 

^3 

1 

2 


■|(u 

+ 0)^ 




u® + 

■ .40 

J 

— J 



(2.60a) 


(2.60b) 
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.(0) u2 _ c2 


2C 

t 2 


[(acf + i)c^ + 

8crC L. -J 




(1) ^ oj^ _ pj2 


80“= 

4^^ = (u + C)2 

* 8C2 

:(1) = (u2 ^ c2) 

4C^ 


4^' = lie 

(U - C) |u2(5u + c) + oc 2(7U - 
X^^^ = ■*■ [^(5U - C) + crc2(7U + c"^! 

Xp^ = ■-- - - g+1 )C^ (U+C ) ^ + U2jc2(4a-3)(a+l) + 6UC(a+l) + u2(Sa-^^ 

= |(<?+l)C^ - 8UC + 7(a+l)U^ + 

r— ^ 


> 


x(2) 

^5 


U^j(5a-3)U^ + 4(a+l)UC - (a+l)C^'| 

x(2) ^ ^ ^ 

32C^ L. 

2 

x!^2) = _ 2 —^ (u + C)(3U + C) 

' 32C^ 

^ ~5^ r(«T + l)C^ + (CT - 7)U^ + 2 (ct + l)uc] 
320 L> 

X^2) = (U - C)|(3a + l)C - (a + 3)^ 

^10 


£ 

64C 


^ (U + C) (3ff + l)C + (a + 3)^ 


X 


(2) _ + l) ^q 2 _ y2j 


11 


16C' 


(2.60c) 
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T 


(O) _ U + C 


T, 


( 0 ) _ 


2C 

U - C 
2C 


m(l) 

■^1 




(1) . 2 ^ (U . c) 


T. 


m(l) 

■^3 


8C' 
g + 1 

8 c 2 


(U + C) 


T 


(l) ^ g - 1 _U 
4 c2 


4 

( 2 ) 


+ ctc2)(2u 2 + aC^) 

^ 48g^C* 

(2) _ 

2 

( 2 ) . 


g 

64gC 


T. 


£ 

64gC3 


T 


T, 


( 2 ) _ 

4 

( 2 ) __ 


.(2) _ 


j^^CC - 3U) + U^(3C - 5U^ 

^ j^^(C + 3U) + u2(3C + 5U^ 

— ^ |gc2|c(a-3) - U(g+1^ - u2j^(7g+3)C + (Sg-3)^ 
^ |^^[(o-3)C + 3(g+l)T^ + U2|j[5g-3)U - (Sg+l)^ 
^ |3(g - l)U + (g + 3)q 

;2C^ L J 


> (2.60d) 


32gC 

a 


,( 2 ) ^ _ 


32C 

g" - 1 

32C^ 


(3U + C) 


|jg - 7)U + 2(g + 1)^ 

j(2) ^ _ (g t _ (,) 


64C 


T. 


T, 


( 2 ) = 

10 

( 2 ) 

11 


_ (g + l)(g + 5)_ ^ 

i3 


64C' 

(q + 1)^ Tj 

3 ^ ‘ 

16C 
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It is to be noted that in equations (2.58) the terms independent of the 
disturbance function are exactly those given by equations (2.45) for the 
steady-flow solution. Thus^ the solution for the unsteady flow also 
generates the underlying steady duct flow below the "rest” characteristic 


To illustrate the character of the solution, an example was computed 
for a duct with an initial exit Mach number of 0.35 (i.e., on x = O) . 

For this Mach number and for a = 5 (t = 7/5) with the flow coining from 
the right, we have U = -0.37879 and C = 1.08226. For the disturbance, 
a "finite ramp" function 



for 


0-$ t < 8 
t > 8 
t < 0 


(2.61) 


was chosen with B = 0.0086905 and 8 *= 4.0. For an € of 0.1, this 
represents about a l/3 percent disturbance in sonic speed on the boundary 
X = 0 at the conclusion of the "ramp." Equations (2.58) were programmed 
for computation by an IBM type- 650 digital computer by using the Bell 
Interpretive System (ref. 9). The program was set up so that the varia- 
tions of the flow variables on lines of constant x could be calculated. 
The values of the coefficients of equations (2.58) (cf. eqs. (2.60)) for 
the initial values given above are listed in table (2.1). The results of 
the computations are shown in figures (2.18) and (2.19) for x = 0.5 and 
1.0, respectively. In both figures the region between a = 0 and p - 0 
is seen to be an unsteady flow field, as expected from the nature of the 
initial conditions. In this region both u and c vary. In the region 
bounded by a = 4 and p = 0, the particle velocity is essentially coS^ 
stant; this reflects the constant value of u between these characteris- 
tics on X = 0. In this region the rate of change of sonic speed" '^th 
time increases above that below p = 0. In the region bounded by a 4 
and p = 4, the magnitude of the particle velocity decreases, and the 
sonic speed continues to- increase until the P = 4 characteristic is 
reached, at which time a new steady state is achieved. 


In both figures the first- and second-order solutions are quite 
similar in shape. The principal difference is in the level of the solu- 
tion which, however, reflects principally the accuracy of the representa- 
tion of the underlying steady flow. The changes in level in going from 
the first- to the second-order solution is greater at x = 1 than at 
X = 0.5, as expected. It is obvious also that the time interval during 
which the particle velocity is almost constant is shorter at x = 1 than 
at X = 0.5, which is to be e3q>ected from the triangular nature of the 
region between P = 0 and a = 4 (shown in the inset in fig. (2.18)). 
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In order to conrpare the perturbation- series solution vith an exact 

solved by the method of characteristics with 
a finite difference technique. Consider a honcharactetlstic arc 1 and 
two points on this arc, J and (j+l) (see sketch). At these 



X 


points the dependent variables u and c are known, io determine the 
location of point (l+l,j) and the values of flow variables at this pointj 
equations (2.36) are -written as difference equations. The resulting 
equations are: 


"1,J °i,J' 


*1+1, J - -i,J 

^1+1, J - ^l,j+l = ("^i,J+l - ®i,j+l)(%+l,j - -tl.J+l). 

U. 4C. 


("1+1,J - Ui,j) + ^(°i+l,J - Ci,j) = € ^ 

1^ 


•l.ri,.1 

+ c 




1^ tJ 


(x 




- X 




), 


(u. 


1+1, J 


'^l,j+l^ " '"^'"l+l,J - °l,j+l^ 


- € 


^,j+l®i..1+l f . 

^i,j+l ■ °l,j+l “ ^^<5+1 • 


(2.62) 


f ^ solved first; this yields the coordinates 

Of the point (1+1^ j). Using these values^ the last two .equations yield 
the values of flow variables at this point. This is a first approximation 
4 accuracy results when average -ralues, rather than 

initial values, are used for the slopes of the characteristics appearing 
in the difference equations . For this purpose, average values of the flow 
parameters that appear In the expressions for these slopes are computed 
by using the values of j and j obtained from the preceding 
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approximation . These average values are designated by barred symbols and 
are given by; 


^i 

°i 


%J+1 “ I J+1 " ) 

°i,J+l “ 2 [°i,j+l ■*■ °i+l,^ 


(2.63a) 


(2.63b) 


Substituting these values into their proper positions in equations (2.62) 
yields the following equations for the second approximation of the loca- 
tion of point (i+1, j) and the values of the flow variables at this point: 

^i+l,J " *i,J " 

^i+l,j " *^i,J+l “ ^i,J+l ■ °l,J+l^^‘‘^i+l,J “ ’‘^l.J+l^" 

K.1,J - -i,j) ^ ^ - *1,J> 


(u. 




1+1, j “i,J+l 


) - a(c 


i+l,j ■ °1,J+1 


) = 


- € 


-p 7 ^^1+1, J ^i,j+i^- 

^,J+l - °i,j+l^ 


(2.64) 

This iteration procedure is repeated until the desired accuracy is ob- 
tained^ that is, until the difference between the values of the flow varia- 
bles for successive iterations is sufficiently small. 
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These equations were also programmed for solution on the digital 
computer. The characteristic net for the problem is. shown in figure 
( 2 . 20 ). The variations of the flow variables at constant values of x 
were obtained by linear interpolation between the values at the net points 
and are shown as the data symbols on figu^s (2.18) jnd (2.19). In both 
the figures the method-of-cheiracteristics Solution i^ in excellent agree- 
ment with the second-order perturbation-series solution. The differences 
are quite small and are slightly greater at x = 1. .than at x = 0.5. 

Prom this, the perturbation- series solution is concluded to give an accu- 
rate representation of the unsteady flow. 

To compute the series solution (both first- and second-order) at a 
given value of x required about 2.5 minutes of machine time. The solu- 
tion of this same problem by the method of characteristics, using the same 
computing machine, took 12 times as long in machine time alone. The ratio 
of machine times will increase as the magnitude of x Increases, because 
of the hature of the method-of- characteristics solution. To the machine 
time must be added that required to plot the characteristic net and to 
perform the necessary interpolations. The advantage of the perturbation 
solution from the computational point of view is quite marked. Further 
comparisons of this nature will be made in section IV. 

3. The "simple- wave- type" disturbance' . - In the_precedlng paragraphs 
the solution for an unsteady flow was obtained with both dependent varia- 
bles arbitrarily specified oh a time-like arc. As ejipected, this led to 
a solution with an unsteady flow below what normally might be called the 
"rest" characteristic. As noted earlier in this section, in order to ob- 
tain solutions with steady flow below the "rest" characteristic, only one 
datum can be arbitrarily specified on a time-like ardi. That is, there is 
a relation between the dependent variables on such a boundary. It was 
attempted to determine whether a general u,c-relatlon on such time-like 
ar.cs could be derived which would permit specification of these variables 
on the boundary and still produce solutions with a steady duct flow below 
the "rest" characteristic. Unfortunately no such reliction could be es- 
tablished analytically. It was therefore decided to investigate some 
simple u, c-relations to see whether the desired result could be approxi- 
mated. Taking a clue from the nature of the zeroth- order equations, which 
are the same as those for one- dimensional flow, it was decided to investi- 
gate a disturbance in which the Riemann variable Q = u - ac remained 
constant on the boundary. One- dimensional flows in which a Riemann varia- 
ble is constant throughout the flow field are called simple waves. Such 
flows are characterized by the fact that they exist adjacent to a region 
of constant state. Because of the nature of the u,c-relatlon, such 
disturbances shall herein be called "simple- wave- type" disturbances. 

For such "simple-wave- type" disturbances in an "epsilon duct" the 
problem may be stated as follows : Given a steady duct flow with flow 

variables U and C at a given cross section. Let there be a disturbance 
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in the flow at this section, commencing at a specified time, of such nature 
that the parameter Q = u - ac is constant on the boundary. Determine 
the ensuing flow in the duct. 


The boundary conditions for the problem are then (choosing x = 0 
as the cross section and t = 0 as the beginning of the disturbance for 
convenience), on x = 0, 


u =s U + ea f(t),] 
c = C + e f(t). 


(2.65) 


where 


f(t) =0 for 

+ ac2 = or + 1. 



( 2 . 66 ) 


The boundary conditions are to be satisfied as follows; The line 
X = 0 is to be mapped on a = 6 such that a = p = t, that is, according 
to equations (2.43) and - — ^ 

u(0)(a,a) = uC0)o,p) = U, 

= 0 f(a), = cf f(P)» 

u^^^(a,a) = u^^^(p,p) “0, n > 1} 

. . . . > (2.67) 

c^0^(a,a) = c^0^O,p) » C, 

o^^^Ca.a) . f(a), - fO), 

c^“^(a,a) = c^“Hp,P) =0, n > 1. 

J 

For these boundary conditions we see that on 


Q = u- crc = U- aC, 


X = 0 we always have 
( 2 . 68 ) 
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(2.69a) 


> (2.69b) 


that Is, a constant value of one of the Eiemann variables on the boundary. 
For these boundary conditions the solutions up to second order are: 

= U, 

u(l) = - P) + u|^^f(p), 

u(2) ^ - p)^ +-U^^^(a - p)f(p) + f(T)dTj 

= c, - - 

c(l) ^ - p) + f(p), 

c(2) ^ cp^(a - p)^ + - p)f(p) + f(T)dTj 

Hi 

= X^°^(a - p), 

= X^^^(a - p)^ + x|^^(a - p)f(p) + x|^^ P^ f(T)dT, 

J.2) ^ x^^^(a - p)^ + X^^^Cci - P)^f(p) +-X^^\a - p)[f(p^^ + 

P (t - p)f(T)dT + x^^^f(p) P f(T)dT + 

jf(T^^dT + x!^^^ ^ (a - T)f(T)dTJ 


/(2.69c) 




RQ?Q 
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^.(2) = 


- $f + 4^'(“ - P)fCS) + 4^^ r 

. f,)5 + lC2)(a - ^)2f(^) + 4=)(a - P)[f0)]' 

B P 

,j«r (t - p)f(Tr)dT + r f (T)dT + 

(a - r)f(T)dT. 




> (2.69d) 




The coefficients appearing In eq.uations (2.69) are functions of 
C; and o and are listed belov; 


uCi) 

1 


uc 

2 


u. 




u. 


(1) _ 


= a 


(2) _ _ (U - C)(U + oC) 


u, 


1,(2) 


4C 

(U + C)(U + oC) 




4C 


(2.70a) 




(1) Ij2 


2a 


.(2) ^ . ufc 


,( 1 ) 


= 1 


1 

,C2) 


4a 


(U - C)^ 
4C 


2 

tC2) _ CU + C) 


> 


4C 


(2.70b) 
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,(o) 


r(l) 


2C 


- - ^( 5 / [( 2 ° + ^ “!] 


.(1) 


a + 1 

4 c 2 


(U - O' 


M = (u + C )2 


.( 2 ) _ 


r( 2 ) = 


r( 2 ) _ 


4 c 2 
Xj 2 

■ 480C 


|(30 + 2)C^ ^ 

^ |oG^(C - 7U) + U^[j4tT - l)C - 
(g + . A)(5g + , ll (u . C) + 


£ 

32 oC 2 


+ OT. 


( 2 ) 


X 


:p) = .2_±_i . ± 9J . (tj2 + 2aW + crC^) + 

4 16a n2 ' '4 


.( 2 ) _ 


16a 

a2 - 1 


r(2) _ 


8C 

a - 


'■6 


^( 2 ) 


16C 

= |aC^jc + (a-l)^ + j(a-l)U + (a^-a-l)^^ + UT^ 


^ (u + c) + 

- jla - 3)U + (3a - l)^ + 


UT, 


( 2 ) 




X2.70c) 


( 2 j 
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T. 


(O) U + C 




2C 


a?. 


T 


(O) C - U 
2 

Cl) 


2C 


= - iL (u^ + ac^) 


8o^ q2 


T. 


(1) _ cr + 1 

2 “ ^ ■ 

(1) _ a . 1 


H J- (u - C) 
4C^ 


AC 

+ 0C2)(2u 2 + oc2j 

48a^C^ 

^ _ . g^l j^c^(C - 3U) + U^(3C - 5U^ 


(U + C) 


,( 2 ) 


= _ C^ •*• ^)C? _ c) 


16C' 


j(2) ^ U_+C ^ ^ ^ u2(3o - l31 

16aC* L J 


,( 2 ) 


- 1 
80^ 


(U + C) 


A Per - 3)U + (3a - l)3 
16C L. ^ 

j(2) .2^, IHJ:^ (o2^rt32). 


(2.70d) 


To Illustrate the nature of the solution for the " simple- vave- type" 
disturbance, an example vas computed for the same values of U, C, €, and 
or used in the preceding example. The disturbance function f(t) is 
again defined by equations (2.6l) with B = 0.0086905 and S = 4.0. The 
values of the coefficients (eqs. (2.70)) for these values of the parame- 
ters are given in table (2.2). 

The results of the IBM computations are shown in figures (2. 21) to 
(2.23) for X values of 0.5^ 1.0, and 2.0, respectively. At aHI values 
of X the character of the variations of both u and c are similar, 
consisting of almost linear changes in the variables between their initial 
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and final values. The most interesting aspect of the xesults is the faclj 
that in the region between a = 0 and ^ = 0 the variables are essen- 
tially constant at their original values. This indicates that the 
"sln^ple-WBve-type" disturbance yields a flow in which the p = 0 charac- 
teristic approximates a "rest" characteristic. The principal difference 
between the first- and second-order solutions is again in the level of the 
variables; this reflects the convergence of the- terms representing the 
steady flow to the exact solution. It is also seen th^at the slopes of the 
second-order solution during the period in which the flow quantities are 
changing are greater than those of the corresponding first-order solution. 
The maximum departure from constant conditions in the region between 
a = 0 and P = 0 occurred in the second-order solution for c at 
X = 2.0. The error in c in this region Ib, however, only 2 percent of 
the total change in the variable during the transient, a quite reasonable 
approximation to steady conditions. Thus, even the second-order solution 
yields an essentially steady duct flow below p = 0. 

Two other points should be noted about the nature of the transient. 
First, as the conpresslve disturbance is propagated upstream, its magni- 
tude increases. For example, at x = 0.5 the total change in u is 
0.0186 or 4.6 percent of its initial value. At x = 2.0 the change in u 
is 0.0235 or about 4.8 percent of its initial value. Corresponding 
changes in c are 0.00366 (0.34 percent) at x = 0.5 and 0.00426 (0.38 
percent) at x = 2.0. The second point of interest is that at each cross 
section the value of Q remains essentially" constant at its original 
value during the transient-; At x = 0.5 the initial value- of Q is 
-5.8042, and its terminal value is -5.8039. At x = 2 the corresponding 
values are -5.8513 and -5.8493. This fact will play an important role 
in the next section. “ 

To check the accuracy of the perturbation-series solution, this 
problem was also solved by the method of characteristics. The character- 
istic net is shown in figure (2.24). Again, the values of the flow varia- 
bles at fixed cross sections were obtained by linear interpolations between 
the net points and are shown as the data symbols on figures (2. 21) to 
(2.23). At x= 0.5 and 1.0, the method-of- characteristics solution is 
in excellent agreement with the second-order perturbation-series solution. 
At X = 2, the agreement is not as good but is still quite acceptable; 
the maximum deviation occurs in the variable c, being about 6 percent 
of the change in c at this station. It is of interest to note that- at 
X = 2.0 the method-of-characteristics solution exhibits a slight decrease 
in c in the region between a =s 0 and p = 0, as does the series solu- 
tion. Thus, over an area ratio of 1.22 the series solution is seen to be 
quite accurate for the unsteady flows investigated. 
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F. Suimnary 

In this section the method of coordinate perturbation was applied 
to the solution of quasi-one-dimensional flov problems. It was demon- 
strated that, for ducts for which the logarithmic derivative of area varia- 
tion with respect to the space coordinate is a function of the smallness 
parameter €, perturbation- series- type solutions for both steady and un- 
steady flows can be obtained. 

> By means of numerical exairples it was demonstrated that the solutions 

j for the steady flow converge to the exact duct flow solution. !Ihe princi- 

pal factors affecting the rapidity of the convergence were shown to be 
the area ratio and the magnitudes of the initial data. In general, good 
representation of the flow was obtained over area ratios as large as 1.8:1. 

For unsteady flow it was shown that solutions could be written for 
arbitrary disturbance functions in both flow variables specified on time- 
like arcs. For such initial data the problems were termed "inverse" prob- 
• lems because, from a physical viewpoint, the solutions give the unsteady 

flow in other regions of the duct (at earlier times) which would produce 
the specified variation of the flow quantities on the boundary. It was 
also shown that a reasonable approximation to the "direct" problem (i.e., 

- one in which a region of steady flow bounded by a "rest" characteristic 

exists) is obtained by specifying that one of the Riemann variables is 
constant on the time- like boundary arc. This disturbance has been termed 
a "simple- wave- type" disturbance. A numerical example indicated that, 
even for the second-order solution, the maximum change of the depend ent ~ 
variables below the "rest" characteristic was of the order of 2 percent 
of the magnitude of the excursion of the variable during the transient. 

For unsteady flow the solutions obtained contained terms involving the in- ' 
dependent variables a and ^, the disturbance function and definite 
integrals thereof. Therefore, solutions for integrable disturbance func- 
tions are easily obtained. 

Comparison of the unsteady- flow solutions with those obtained by 
the method of characteristics indicated that the second-order series solu- 
tion accurately represents the flow. As was the case for steady flow, the 
accuracy of the series solution decreases with increasing distance (i.e., 
area ratio) from the boundary. The series solutions offer an advantage 
in computational effort over the method of characteristics calculations. 
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III. UTOEEIACTION OP A SHOCK WAVE AND A 
DISTURBANCE IN A DUCT PLOW 

In section II^ perturbation- series solutions for unsteady quasi-one- 
dimensional flow were developed. These solutions are now applied to the 
determination of the motion of a shock -wave in a diffuser flow under the 
influence of a disturbance originating downstream of the shock. An ap- 
proximate relation is used for the shock transition. It 'is shown that, 
for a linear variation of the flow variables in a "simple- wave- type” 
disturbance, an analytical solution for the shock path can be obtained to 
first order in € . This solution is compared with both first- and second- 
order solutions obtained by numerical integration of the differential 
equation of the shock path. 


A. The Equation of the Shock Path 

Por weak and moderately strong shock waves, the_ shock velocity can 
be approximated by 

V= (^s+^s) +i [(^s+^s) - (^s+=s)] [(’^s+^s) - (^s+^e)]^- C^-l) 

0 

where V= v/c^ (cf. ref. 6, eq. (72.06). In this equation the subscript 
s refers to the immediate shock location, the barred quantities represent 
conditions upstream of the shock, and the unbarred terms represent condi- 
tions downstream of the shock. This equation is correct to second order 
in the shock strength. Within the same degree of accuracy, the shock 
transition is homentroplc; and, for a "forward- facing" shock wave (when 
the fluid particles enter the shock from the side of larger x- values), 
the Rlemann variable Q is constant across the shock; that is, 

Ug - (TCg « Ug - crcg. (3.2) 

These equations are applicable for weak ani moderately strong shocks^ and 
the development that follows is valid only for such shocks. 
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Consider a diffuser flow in an "epsilon duct" -witli a stationary shock 
wave at the cross section x = X (see sketch) . 



For a diffuser flow with an Intlally stationary shock, the conditions 
upstream of the shock (i.e., u and c) are functions of x alone and 
are not influenced hy any subsequent shock motion (unless, of course, the 
shock is regurgitated). This, and the approximate shock relations (eqs. 
( 3 . 1 ) and ( 3 . 2 )), imply that, regardless of the shock velocity (within 
the limits of accuracy of eq. (3.1), of course), the value of Q remains 
constant at a given duct cross section as the shock passes the section^ 
although Q changes from section to section. As noted in section II- E, 
a "simple-wave- type" disturbance has the property that the value of Q 
at a given section remains approximately constant at its Initial value 
during the transient. Thus, if the disturbance downstream of the shock 
is of this type, it will approximately fulfill the boundary condition at 
the shock, that is, equation (3.2). This is equivalent to neglecting the 
"back reaction" of the shock on the flow field behind the shock. In other 
words, the flow field behind the shock is completely specified by the 
solution for the "simple- wave- type" unsteady duct flow. The subsequent 
development is, therefore, based on the "simple- wave- type" disturbance. 


The velocity of the shock wave is given by 

dx/dt = V. 


(3.3) 


If the flow field behind the shock is described in terms of the character- 
istic parameters a and 3, we have 


Thus equation (3.3) 


dx = x^ da + Xp 
dt = tj^ da + tp 

can be written 



3 ^ da + Xp d|3 
t^_ da + tp dp 


(3.4) 


(3.5) 
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The partial derivatives of x and t are related by the equations of 
the characteristics; that is, 


= (u * 0)^1 

Xp = (u - c)tp. J 

Substituting equations (3.6) into equation (3.5) and solving for da/d 3 , 
we obtain: 


§ 2 : = - ^^s - °s)3^b 

CCUg + Cg) - V]t"^ 


(3.7) 


■where the subscript s again indicates that u and c are to be evalu- 
ated immediately downstream of the shock. Substitution of the expression 
for V (eq.. (3.1)) then yields 


da 

[2 ^^8 ~ '^S^ 2 ^ 3 — 

[(Ug + Cg) - (Ug + Cg) 


d 8 " j 

[2 '^S^ “ 2 ^’^s °s^ " 

3 

[(Ug + Cg) - (Ug + Cg^ 

2 ) 

i 

1 a 


(3.8) 

In this equation, t and the "immediately downstream" quantities Ug, Cg 

are known functions of a and 8 given by equations (2.69) for the 
"simple- wave- type" disturbance. In the flow field behind the shock, up 
to and including the instantaneous position of the shock, x is also a 
known function of a and 8 . Therefore_^ in principle at least, the 
"immediately upstream" values Ug and Cg can be considered as functions 

of the characteristic variables as u and c are functions of x alone. 
Thus equation (3.8) may be written as 

^=G(a, 8 ). (3.9) 

This is an ordinary differential equation for the relation between a and 
8 on the shock path. Substituting the solution of this equation in the 
equations for x and t (eqs. (2.69c) and (2.69d)) yields a parametric 
representation of the shock motion resulting from the disturbance specified 
at the boundary. 

Equation (3.8) can, of course, be integrated numerically for specific 
cases. It is of Interest to see if an analytical solution can be obtained, 
and such a development follows. 
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B. Approximate Solution for the Shock Path 

The first-order solution for the unsteady flow field behind the shock 
for a "simple- wave- type” disturbance is repeated below for convenience: 


u = U + eju^^^Ca - ^) + 
c = C + e|c^^\a - p) + 

X = - p) + ejx^^^Ca - p)^ + - p)f(p) + 

f (T)d^ , 

t = Tp^P + - p)^ + - P)f(P) + 

f(T)d^ . 


(3.10) 


From the last relation it is easily shown that 

= ^2°^ + ejaT^^^Ca - p) + T^^^fCp) 

^P = + ejjsT^^^Ca - p) + - 

- P)f(p^. 


T^^hf(P) + 



(3.11) 


The first step in the development is to express u and c in terms 
of the characteristic parameters. For a steady duct flow with conditions 
U and C specified at some area Aq^ the following relations hold: 



+ ac^ = (a + l) . 


(3.12) 


(3.13) 
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Assume that at some area A(x) = A, upstream of and smaller than Aq, but 

still greater than that at -which sonic flow would exist, a shock occurs 
(see sketch on p. 47). At this station we thus have on the downstream 
side of the shock 


Aq 


a/2.„a 
a ^ UC , 


u(X){(o + 1) - [^u(X^ 


vT12| 




For a stationary shock the conditions on its upstream side can be deter- 
mined frcxn equations (3.1) and (3.2) with V = 0, or.jnore exactly from 


UgUs = 1, 


Prandtl equation; 


(3.14) 


Ug + uCg = Ug + ac^, energy equation. 

The steady supersonic flow ahead of the shock is now given by 


( 3 '. 15) 


A _ • u (X)[je(X)_ 


A 

u 


[(a + 1 ) - u^ 




(3.16) 


+ ac^ = (a + 1 ). (3.17) 

Thus^ specifying an initial steady shock location by x = X and the 
conditions U and C at the end of the diffuser^ Aq, completely deter- 
mines the variables u and c ahead of the shock. 


In the problem being considered^ the duct-area variation is again 
•written 


a/a = 


-e(x-X) 

e 


Equation (3.16) therefore may he written 



(3.18) 


(3.19) 


where, for simplicity of notation, u(X) and c(X) immediately on the up- 
stream side of the shock at its initial location have been replaced by 
U and C, respectively. 
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In order to develop the upstream flov variables as functions of a 
and P, it is assumed that these quantities can be represented to first 
order by 


u = U(l + 
c = C(l + 



(3.20) 


where and axe functions of x. Expanding equation (3.19) to 

the first power in e after substituting equations (3.20) yields; 


J?72 

OC + € 


2( . yj-) . (1 - = 0C^|1 + e I (3 



Upon substituting for x from equation (3.10), the right hand side 
becomes 

00^|l + e f [ca - ^) + 0(e) - A, 

whence 



cr 17 

“ + ^ (1 - U=) 



(3.21) 


Substituting equations (3.20) into (3.17), expanding to first order, and 
using ( 3 . 21 ) gives 



1 if 

1 (1 - U2) 



P) 



( 3 . 22 ) 


A term that appears in the shock velocity relation is (u + c) . This can 
now be written in terms of a and ^ as 

uH.;.(o^c)^e 'B5(n-°c) (3.33) 

(cr + 1)(1 - U^) L 2C J 

Thus the flow variables ahead of the shock are given to first order as 
functions of the characteristic variables. 
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We now turn to the evaluation of the right hand side of equation 
( 3 . 8 ) to first order. Consider the bracketed factor in the numerator j 
that is, 


R = i ( 3 C 3 - U 3 ) + I (Ug + Cg) + ^[(^S + °s) - (""s + °s^^‘ (3-24) 

Substituting the expression for the flow variables from equations (3.10) 
and ( 3 . 23 ) and expanding equation (3.24) to first order in € gives 

R = Rq + e jR^(a - S) + RgfCp) + R 3 (3.25) 

where - 


R^ = 


R, = 


3C - U . U + C 

^ ^ T 




+ c) - (u + c) 


3^^ 


_ U(5U+aC) ^ UC(U-gC)(U^-c2) ^ [ (u+C) - (u + C)3u(aC-U) _ 
4 (ct+1)(1-U^)C ' ~ 


4a 


8oC 


[(U+C) - (U+C)]U(U-gC)Cu2-c2) _ [(U+C) - (U+C)]2 U^(u2-c2) 
8(a+l)(l-u2)c 16(a+l)CC(l-u2) 


1 = 3 - P ^ (g + 1)[(U +_C) - (U + C)] 

‘2 2 4C ' 

. ^ UC(U-t^) [(U^) - (U+C)3u(U-gC) _ [(U+C) - (U+C)]^u2 

^ 2(g+l)(l-U^) 4(g+l)(l-U^) 8(g+l)(l-U^)C j 


(3.26) 

Similarly, the bracketed factor in the denominator of equation (3.8), 

S = I (Ug + Cg) - I (Ug + Cg) - [jUg + Cg) > (^g + ^ ^ ^ 

reduces to -■ ■ — 

S = Sq + €js 3 _(a - ^) + SgfXp) + S 3 (3.28) 
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vhere 


„ U + C U+C [ (U + C) - Cu + C)]^ 

^ 2 ^ 

g _ U(gC-U) UC(U-oC)(U^-C^) [(U+C) - (u+C)]u(aC-U) ^ 

■ ^ 4(0+1) (l-u 2 )C 




8aC 


[(U4C) - (U+C)3u(U-aC)(U^-c2) ^ [(U+C) - (U4<3) ]^(U^-C^) 

8 ( 0 + 1 ) (l-U^)C 16(0+l)C(l-l^)C 

cj g + 1 (0 + i)[(u + c) - (ir + c)3 

2^2: 4C 

g ~UC(U-0C) ^ [(U+C) - (U+C) 3 U(U-oC) ^ [(U+C) - (U+C)] 2 u 2 

^ 2(0+1) (1-u 2) 4(0+1) (l-ij 2 ) 8(0+1)C(1-u2) ^ 

(3.29) 

Substituting these expressions and those for the psurtiaJ. derivatives of 
t into equation (3.8) yields --- 


to _ {Bq + e[R3^(a>-p) 1- Bgf(^) -I- 

{S(J + €[Si(o-^) + SgfCp) + + el 2 T(l)(o-p) + - 4 ^^f(a)I- 

which, when expanded to first order, becomes 


(3.S0) 


g: = Bq + e|j:a-p)B3_ + f(3)B2 + (a-p)f‘(3)B3 + f(a)B^ + XB^ j (3.3l) 


wliere 




Bq = BoTi yso4°^ 


,C0) 


Bi = 

SqT^ ^ L -* 

®2 = TZm ■" ^?^®2 - + 4° ^Sg)], 

®5 = Vi^V So4°’> 


/ (3.32) 


=4 = V' 


S3 

~ ^0* 


J 
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Equation (3.31) cannot be formally Integrated unless the disturbance 
function f(r) is known. If f(T) is a linear function of its airgument, 
equation (3.3l) becomes a linear differential equation and can be easily 
integrated. Making use of this possibility, we introduce for the disturb- 
ance, the function 



(3.33) 


Because of the nature of this function there are three zones in the 
x,t-plane in which different versions of the right hand side of equation 
(3.3l) hold. These zones are bounded by the characteristics shown in 
the following sketch; 



In the three zones, the function, its derivative, and integral assume the 
following values: 


Zone 

f(a) 

f(P) 

f'(p) 

j f (x)dT 

I 

Ka 

K^ 

K 

1 K(p2 - a^) 

II 

K8 

KP 

K 

Kjj p2 + 1 g2 _ 

III 

K8 

K6 

0 

KB(p - a) 
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o> 

to 

CVJ 




In zone I, equation (3.3l) becomes 


da 

dp 



Bq + 




(3.34) 


The solution of this equation is 

eDgP 

^ ^1^ ^^s ” ^0^^ ^ (3.35) 

vhere 


(Bq + eXBg)! 


+ KCBg + B^)] 

+ 

|k{E3 - B3) - b‘ 

6 

B3_ + K(B3 + B 4 ^ 

2 


K(B2 - B3) . B^^ 

+ K(B3 + 

Dg = B^_ + kCBj + B^^); 

and use has been made of the initial condition that, -when p = 0, a *= . 

The value of is determined from the initial shock location. Relations 

of form similar to equation (3.35) hold for the other zones; that is, a 
is the sum of a linear function of P and an e 2 q)onential function of p. 
Substitution of equation (3.35) into the expressions for x and t in 
equation (3.10) yields, to first order in e, a parametric representation 
of the shock path under the action of a "simple-vave-type” disturbance of 
the form given by equation (3.33). The shock path cannot be given ex- 
plicitly in X and t because the nature of the a, P-relationship along 
the path does not permit the elimination of a and P algebraically. 


> (3.36) 




C • Examples 

The solution for the shock path presented above is correct to first 
order in €. As noted previously, the expressions for the shock transi- 
tion used in deriving this result are correct only for veak or moderately 
strong shock -waves. Because of this fact, shocks occurring at an upstream 
Mach nuiDber of about 1.5 (for y = 7/5) are about as strong a shock as 
can be accurately represented by the approximate equations. Therefore,, 
the examples chosen will have an initial upstream Mach number of 1.5. 

Another fact that enters into the determination of the region of 
applicability of the approximate solution for the shock path is the 
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accuracy of the representation of the initial steady flov offered by the 
approximate relations used. The flow conditions upstream of the shock 
wave have been represented by equations (3.20). For a shock at Mach 1.5 
we have U = -1.36458 and 7 = 0.90972. Using equations (3.2l) and (3.22) 

to evaluate and the variation of^- u and “c with (x - X) has 

been computed for € = 0.1 and is shown in figure (3.1), along with the 
exact values of these parameters. For values of (x - X) up to 0.4, the 
firstj^order representation of the flow variables is quite good. At 
(x - X) = 0.408, the Mach number of the flow is 1.425. Therefore, ii^the 
disturbance magnitude is limited to that which would result in a finaJL 
steady- state shock location at an upstream Mach number of about^ 1.4, the 
representation of the upstream conditions should be adequate. 

The disturbance function is specified on an x « constant line down- 
stream of the initial shock location. From the examples given in section 
II for the steady-flow solution, it was evident that the first-order rep- 
resentations of the flow variables are accurate only for relatively short 
distances upstream from the boundary. Because of this, the farther from 
the initial shock location the place x = constant is selected, the 
greater will be the error in the shock path computed by the first-order 
solution. To illustrate the accuracy afforded by the first-order solution 
for the downstream flow variables, two examples have been computed. For 
an upstream Mach number of 1.5 the Mach number on the downstream side of 
the shock is 0.7011. If the boundary on which the disturbance is speci- 
fied is located at the section at” which the initial Mach number is 0.7, 
we have X = 0.007, U = -0.73179, and C = 1.04541. For these conditions 
the downstream flow variables are shown in figure (3.2). (it will be noted 
that X does not enter this calculation except" by determining U and C 
at X = 0.) For values of x up to 0.4 the first-order representation 
is quite good. If we choose the boundary at the section where the initial 
Mach number is 0.6, we have X = 0.830, U = -0.63481, and C = 1.05801. 

This case is shown in figure (3.3). At the initial shock location, the 
first-order values of the flow variables are still reasonably accurate. 
Upstream of this station the departures become quite large, especially 
for c. It should be noted that the values of the downstream flow varia- 
bles upstream of the initial shock location (which do not exist prior to 
movement of the shock) serve as the basis for the new downstream flow as 
the shock moves upstream. Hence, errors in final shock location will be 
larger when the boundary conditions are given farther from the initial 
shock position . _ 

The first shock path calculated is for a shock initially occurring 
at a Mach number of 1.5 in a duct with €- = 0.1. The disturbance function 
is specified at the downstream section at which the initial Ifech number 
is 0.7 (X = 0.007) as in figure (3.2). The values of U, C, U, and C 
have been given above. The magnitude of the disturbance was selected so 
that the final shock position would be that corresponding to an upstream 
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Mach, number of 1.43. For a 5 of 4, the corresponding value of K (cf . 
eqs. (3.33)) is 0.010025. For this disturbance the shock should come to 
rest at x = 0 . 388 ( cf . the sketch on p . 54, point L) . The coefficients 
of equation (3.3l) are: 

Bq = 1.04087, B 3 = -16.92621, 

B]^ = -0.78888, B 4 = 0.35140, 

Bg = 16.10895, Bg = -1.52893, 


and ttg = 0.02629. 

The equations relating the characteristic variables on the shook 
path are: 

Zone I: a = -1.39236 + 1.172773 + 1.41865e"° *095503^ 

Zone II: a = -1.32790 + 1.168493 + 1.35689e"®'°95®^P, 

Zone III: a = 1.34140 + 3 - 1 . 467780 ***^ 

Substitution of these values into the expressions for x and t in equa- 
tion ( 3 . 10 ) yields the shock path shovn by the dash-dot line in figure 
( 3 , 4 ). The shock asynptotically approaches a final position of x = 0.349^ 
vhich is about 90 percent of the excursion it should have undergone. The 
greater part of the shock motion occurs in zone HI, vhich starts at 
X = 0.076 on the shock path. The value of the final shock location can 
be obtained directly from the a, p-r elation in zone HI. From this equa- 
tion, (a - ^) is seen to have a limit of 1.34140 as P In zone III, 

the integral appearing in the expression for x is a function of (a - p), 
and therefore x is a function of this quantity alone in this zone. 

To check the accuracy of the analytical expression for the shock 
path, eciuation (3.8) vas integrated numerically by using both first- and 
second-order representations of the flov field behind the shock (i.e., 
eqs. (2.69)). The coefficients for these equations are listed in table 
( 3 . 1 ). The numerical integration vas programmed for the digital computer 
by using the Runge-Kutta method of integration and an increment in P of 
2. (Halving the increment of P did not change the results.) The results 
of the numerical integration are shovn as the solid and dashed lines in 
figure ( 3 . 4 ) and are tabulated in table (3.2). For this case the analyti- 
cal approximation congjares very veil vlth that obtained by numerical 
integration . 
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The shock path was then calculated for the same initial shock condi- 
tions vith the disturbance specified at the section at vhich the Mach 
number was 0.6 (X = 0.850 as in fig. (3.5)). For a 5 of 4, the value 
of K for a final shock Mach number of 1.43 is 0.006225. For this case 
the coefficients of equation (3.31) are: 

Bq = 1.33816, Bg = 15.60186, B4 ="0,63237, 

Bj_ = -0.48838, B3 = -15.17776, Bg = -1.56876, 

and ttg = 2.69885. 

The a, ^- relations on the shock path are now; 

Zone I; a = 0.57893 + 1.17443P + 2.11992e“°*®^'7^°P, 

Zone II: a = 0.73649 + 1.166586 + 1.95260e"^-^®296^ 

Zone III: a = 5.08497 + 6 - 2. 5964 7e"0 -04:8846. 

The shock path calculated from these equations is shown in figure (3.5) 
as the dash- dot curve. The curve is similar in shape.. to that of the pre- 
ceding figure and has an asymptote of x = 1.473. The final shock posi- 
tion should be x = 1.211. Thus, specifying the disturbance farther from 
the initial shock location has increased the error and caused the error 
to change sign. Numerical integrations of equation (3.8) for this ease 
yielded the other curves on this figure. (The coefficients are listed 
in table (3.3) and the results of the integration in table (3.4).) For 
this case it is seen that there is an appreciable difference between the 
analytical approximation and the. corresponding first-order numerlcsil 
integration. This difference results from the higher order terms that 
were neglected in obtaining equation (3.31). In the first example, these 
terms were of small import because the disturbance was specified so close 
to the shock. The second-order numerical integration yields a final shock 
position of X = 1.162, which represents a shock movement of 0.332 units. 
The second-order integration in the first example gave a shock movement 
of 0.343 units, which is quite close to the result"-of the present example. 
This reflects the improvement of the representation of the downstream 
flow afforded by the use of the second-order representation of the disturb- 
ance function. 

In these examples the shock displacement obtained from the second- 
order numerical integration (the most accurate result) was always less 
than that which should have been obtained. Thljs error _in shock displace- 
ment results from the fact that the values of U and '”C used in the 
calculations were obtained from the exact shock transition relations (eqs. 
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(3.14) and (3.15)) rather than the approximate relations used in deriving 
the equations for the shock path. In figure (3.6) the Mach nuniber down- 
stream of the shock is shown as a function of the Mach number ahead of 
the shock for both the exact and approximate transition relations. For 
the same downstream Mach nun±>er^ the approximate transition requires a 
greater upstream Mach number than the exact transition relation. Because 
the solution for the shock path employs the approximate shock relations, 
the shock should come to rest at a higher upstream Mach number, that is^ 
after a smaller excursion, than would be indicated from the exact rela- 
tions. If a lower initial shock Mach number were chosen, the error in 
final shock location would be smaller. 


D. Discussion 

The solution obtained for the interaction of a shock wave and a 
” sing) le- wave- type” disturbance in a duct flow yields a shock path that 
may be described as an "exponential decay” approach to a final position. 
This result is quite similar to that obtained by Eantrowitz in reference 
2, wherein he obtained a first-order linear differential equation relating 
X and t on the shock path by linearizing the equations of motion and 
the shock transition relations. 

In the exan^les given above, it was noted that the majority of the 
shock motion occurred in zone Illj that is, where ^ ^ 6. In this region, 
as was required by the boundary conditions, the flow variables downstream 
of the shock have reached their terminal values. Thus, the effect of 
neglecting the "back reaction” of the shock motion on the downstream flow 
conditions is quite pronounced in zone III. This assumption, which is 
inherent in the development presented', is not as restrictive as might ap- 
pear at first glance. In the next section, it will be shown that the con- 
ditions specified on the downstream boundary are reasonable approximations 
for the actual process. 

From the examples presented, it can be concluded that the first-order 
analytical solution for the shock path offers a reasonably accurate rep- 
resentation of the path, provided that the choice of the initial conditions 
are such that the first-order representations of the initial steady flow 
(both upstream and downstream of the shock) are accurate. That is, the 
boundary on which the disturbance is specified should be fairly close to 
the initial location of the shock. Also, the disturbance magnitude should 
be selected so as to result in a shock displacement that is within the 
region of accuracy of the first-order representation of the upstream con- 
ditions. If the boundary conditions are required to be specified at 
fairly large distances downstream of the shock, numerical integration of 
the differential equation should be employed. 
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IV. THE FLOW PIEU> BEHIED.. A MOVING SHOCK WAVE 

In the preceding section a solution for the Interaction of a shock 
wave with a "simple- wave- type" disturbance in a duct:-flow was presented. 
This solution was developed on the basis of two assuiiptions : First, that 

the simple-wave relation (u - cc) = constant held on the downstream 
boundary. Second, that the "back reaction" of the shock motion on the 
flow could be neglected. The solution obtained indicated that under these 
assumptions the greater part of the shock motion occiurred after the down- 
stream flow had reached its terminal state. This leads one to inquire 
how well the results obtained under these assumptions represent the actual 
process. That is, in the flow field behind a moving shock wave, does the 
u,c-relation at a given cross section approximate constant, and is 

the back reaction of the shock sufficiently small in magnitude that it 
may be neglected? 

The problem at hand is, then, to find the flow field behind a shock 
wave moving in a duct- flow. This falls into the category of an "inverse" 
problem as defined in section II and is termed the "inverse- shock" problem. 
The solution to this problem can be viewed as yielding the nature- of the 
disturbance that would result in a prescribed shock path. 


A. Derivation of the Solution 
Consider a shock path in the x,t-plane given by 

X = € f(t) . (4.1) 

The shock velocity is then given by 

V = ^ = 6f'(t). (4.2) 

Because an initially steady flow in the duct is desired, the initial 
velocity of the shock wave must be zero. Therefore, it is prescribed 
that 


f’(Of) = 0 


(4.3) 


and, for convenience, 


f(*b) s= 0, 


t< 0 


(4.4) 
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Aliead of 1)116 shock there is a steady supersonic flow which, for the 
"epsilon duct," can be represented to first order by (cf. section III, 
eq.. (3.20)) 


u 

c 


U(1 + 

C(1 + ec^^^), j 


(4.5) 


vhere U and C are the flow conditions immediately upstream of the 
shock at its initial location^ x = 0 in this case, and 




g 

^ 1 (1 . u^) 

1 „ 


(4.6) 


As shown in section III, these relations yield quite accurate representa- 
tions of the upstream flow conditions for the magnitude of shock displace- 
ment considered. 


After specifying the shock path and the upstream flow field, the 
neict step is to determine the conditions on the downstream side of the 
shock. From the approximate shock transition relations (eqs. (3.1) and 
( 3 . 2 )), the conditions on the downstream side of the shock can be deter- 
mined as follows; 


It is assumed that on the downstream side of the shock the flow 
variables can be represented to first order by 


Ug = U(1 + 

= C(l + ec^^h. 


(4.7) 


where U and C are now the flow conditions immediately downstream of 
the shock at its initial location. Substituting equations (4.7) and 
( 4 . 5 ) into equations (3.1) and (3.2) and expanding the resulting 
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easpression to first order In e gives, from equation (3.1), 

V = I (U + C) + I (U + C) + j^(U + C) - (U + 2 ^ ^ 

6^^ + ^s^^)(3C + U-U-C)- 

(U + C - U - C)^ + 

80 ® 

^ + cc^^bCc - U + C + U)|, ^ 

and, from equation (3.2), 

U - oC + €(Uu^^^ - 0 ^^^^) = U - oC + e(Uu^^^ - oc4^^). 


> (4.8) 


(4.9) 


Substituting equation (4.2) into equation (4.8) and equating the coeffi-r 
cients of like powers of e in equations (4.8) and (4.9) gives 


and 


_ 

U + C + U+ C = --4:CU + C - U- 0)2, 

40 

U - oO = U - oC, 


f ’(t) = i (l3u|l) + ^(^))(30 + U - U - 0) - 

4 (U + C - U - 0)2 + 

4 + CcW)(o . u + u + 0), 


(4.10) 




These equations determine U, 0, and the functions 
these are given by 


(4.11) 


u(l) 


and c 


( 1 ).. 




(4.12a) 
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U = aC + U - aC, 


(4.12b) 


f • (t) X :r 5 — \ aC(U+ C) + 

(cr+l)C(C-U+U+C) (cj+ 1)2(1- u2)c ( 

2C(oC - U)(5C + U - U - C) + UCU + C - U - C)2\ 

2(C - U + U + C) J 


f + X <C(U+C) - 

® (cr + l)U(C- U + U + C) (ct+ 1)2(1- u2)u\ 

2C(gC - U)(^ + U- U-C)+U(U+C-U - C)2 j ^ 

2(C - U + U + C) (‘ 


Substituting equations (4.13) into equations (4.7) gives expressions of 
the following form for the flow quantities on the downstream side of the 
shock: 

= U + eE_f *(t) + eE.x, 1 

> (4.14] 

Cg = C + €E2^f’(t) + sEgX.j 

The are given by: 


E, = 


1 (a + l)(C - U + U + C)" 


Eg = - 


(a + l)2(l - u2) 


orC(U + C) + 


2C(oC - U)(5C + IT - U - C) + U(U + C 
2(C - U + U + C) 


U - C)2^ 


(o + l)(C - U + U + C)' 


Ca + i)^Ci - Jfi) 


^c(ti + c) - 


2C(gC - U)(5C + U - U - C) + U(U + C 
2(C - U + U + C) 


U - C)^ 
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Finally, substituting equation (4.1) into equations (4.14) gives 


Ug = U + eEgf >(t) + e2E^f(t), 


^8 


= C + eE^f'(t) + e%l 2 f(t), 


(4.16) 


These are the conditions on the downstream side of the shock and are, 
therefore, the boundary conditions to be satisfied on this side of the 
shock path. It' is of interest to note that the effects of shock velocity 
appear as terms of first order in €, whereas the effects of~ shock dis- 
placement appear as terms of second order. 


Having obtained the conditions on the downstream side of the shock, 
we now proceed to determine the flow field downstream of the moving shock. 
A solution of equations (2.37) is sought under the boundary conditions 
given by equations (4.16). The boundary conditions are to be satisfied 
in the following manner: The shock path is to be mapped on oo = 6 such 

that t = a = 6; that is. 


x^^^(a,a) 

x^^^(a,a) 

x^^^(a,a) 


x^°^(P,P) = 0} 

f(a), x^^\p,6) = f(p); ? 


x^“^(p,p) = 0, n>ljJ 


(4.17) 


t^°^(a,a) = a, t^°^(p,p) = Pj 
t^“^(a,a) = t^“^(p,p) =0, n > 0; 


(4.18) 


and 


u^®)(a,a) = uC®)(p,p) = Uj 

J^^a,a) = Egf '(a), u^^\p,p) = Ejf '(p)j 

J^^(a,a) = J^\p,p) = E4^f(p)j 

u^“^(a,a) = u^^^(p,p) =0, n > 2j 
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c^®^(a,a) = = Cj 

c^^^(a,a) = >(a), = E^f '(^) 

c^2)(a.^o,) = EgfCa), c^2)(p^pj ^ EgfXp); 

c^^^(o.,a) = c^^^(3,P) =0, n > 2. 

The solutions up to second order are: 




> 


J 


uCo) = u, 

- 3) + '(p) 

u(2) ^ _ pj2 ^ 

U^2)(a - p)f(p); 



cCo) = C, 

c^l) = _ c^l)(a - p) + c|^^f ‘(p), 


> 


cC2) = _ c^2)^^_Pj2 ^ c^^^f(a) + C^^^f(p) + CpUa^P)f ' CP) J J 


t(0) = t(0)p _ T^°^a, 

= - t[^^(o-P) 2 + T^^^jf(p) - fCot^ + T^^^(a-p)f ’(p), 

t|^^ f(T)dT + 

j(2) dT - _ j(2)^ c^p)f(a) - 

T:^^^(a-p)2f(p) + T^2^f‘(p)p(a) - f(p)]j 


t^^^ * “ TpUa-p)^ - _ 


(4.20) 


(4.21a) 


(4.21b) 


> (4.21c) 
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- p). 




,(!) = 
xC2) = 


- x(^Ha-3)^ + x|l^f(P) + X^^^f(a) + x|^)(a-p)f (p), 

- Xp)(a-P)^ + x|2)(a,.p)ji'(p^2 ^ ^^2) f (T)dT 

^^|f'(T^2 dr + xp)(<i.5)i.(p) . x|2)(».p)f(a.) 
x|^^(o-s)^f '(e) + x^^'f '(e)jf(a) - f(ei . 


7 (4.21d) 


J 


The coefficients are given by (use has been made of the fact that 
Eg = 0^E2_) : 


= uc/s 

4^^ = 

y(2) ^ U(aC^ - U^) 
1 8a 


(2) E, n 

4 = ^ + C) + •! Ce^ - CE 2) 

= I + '^^2) - Sc + u)(u + c) 

= Sc 


> 


,( 1 ) 

■^1 






= Es/° 


c(2) „ 


(2) U „ (U -f C)^ i _ V 

9. “ -*^3 ArrP “ 9rt “ ^9* 


“3 ^4^ 
^2 


c(2) _ ™ (U + C)“ U. + i fE 4. nS 1 

^3 “ ^3 ^oC - ^ ^ '.^<1 + 


20 '“4 
20 " 2o '•“4 



(U - C)^ 
4oC 




(4.22a) 


(4.22b) 
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U + C 
2C 

U - C 
2C 


& + u^) 


^ (g + 1)(U - c) 


= (gc2 + u2)(ac2 + 2 u2) 


480^0^ 


2 U + lj(,g + 

^3 1602c3 


(u - c) 


160^03 


jj(0+l)U^ + 0(3-50)UC^ + (40^+50-3)U^C + 0(3-0)C^ 


E ^ - 1 E (g - 1)U ■ 0 - 1 IT 
^4 40C " ^2 4q2 ^ 40 C 


jlo+3)U^ + (20^-0-l)U^C - 0C30+1 )Uc 2 _ cj(o-l)c3j - 

(E, + aEj t F[C<^---Ug±.M 

^ 2 8oc2 8oc2 


EgCU + C) 
1602c3 


jgCg + l)C^ + 4:Ol]C + (30 - l)U^ 


(E, . cE,) uK» , i)c - . m i 

^ ^ 80C2 80C2 


3 2 

32 0^^' 


j^^(C - 3U) + U^(3C - 50r^ 


p2 - 1) IT + C 

■‘3 Qrr2 n3 


P g + 1 
^3 40C2 
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r(0) _ - C2 


20 

t2 


= E 


g - 1 
^ 4crc2 


(U + 0)2 - 


20 




y(l) _ U + 0 „ 

^3 "20 " ® 


'3 

* 4ao2 


= E- 0-+^ (U - 0)2 


r(2) _ _u! 


4a0' 


X 


( 2 ) _ 


48a‘‘ 


|(g02+u2)(ac2+2u2) + a0^|(3a+2)0^ + 


16a2c^ 


(g + l)(U 


C)jT3g 


+ 1)0 - (cj + 


3)^ 


3^^ = j^Qg2g3 |^3a+7 )0 + 5(a+l)l^ - crUC^ |73g-l)0 + (Sg-3)uj 
u2o [(40^ + 3g - 1)0 + (4g2 + 50 - 3)1^ + 


00 ' 


E 


(0 + 1)0 + (3 - 0 )U >+ 


-he 




■4 20 O 


0-1 
'2 4c2 


(U2 + c2) - 


0-1 

40 


u2 


E" 


X 


(0-l)(U+O)f(30-l)O + C0-3)ul - (0-l)(U+O) 

l_ J 4crC^ 


( 2 ) ^ 

1602 c2 'L” " 40o2 

E,U 


rC2) _ -^3 


j^ 02 o^ + (202 + 0 + 1 )UC^ - 2(g^ + l)u2c - ( 0 + 3 ) 1 ^ + 
(®4 + <"^2^ + (g-l)UO + 2U^ 

|0(g+l)O^ + 3g(g+l)UC^ + (50+1 )u2o + (30-l)U^ 4 
Hr (H^J - ^ [(=+1)0^ - C=>-1)DC +-- 2 ;^ 
^7 ^ ” 32o2 ^ **■ -h ^4a-l)C^ - 4UC + 


-.J- - g^) . !■ iUtl (U -_C) 

"^3 


‘8 


=■3 2 

80 ^ 


40 




J 

( 4 . 22 a) 
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These equations give the flow field behind a shock moving on an arbitrary 
path in an "epsilon duct.” The only conditions placed on the shock path 
are those of equations (4.3) and (4.4) and the fact that the shock does 
not leave the supersonic portion of the duct. 


B. Examples 

1. The "exponential” path . - The shock paths obtained in section III 
vere described as "exponential decays" to their final position. There- 
fore, the first shock path investigated was chosen so that it would have 
an "exponential- decay” approach to a final position. Such a path is given 

by 


f(t) = 




t< 0. 


, t ^ 0} 


(4.23) 


This path obviously satisfies equations (4.3) and (4.4). The coefficient 
B determines the magnitude of the excursion of the shock, and b (some- 
times called a "time constant”) determines the rate at which the shock 
approaches its final position. For convenience, the path given by equa- 
tion (4.23) will be referred to as the "exponential” path. 


The first case investigated was for a shock initially at a Mach num- 
ber of 1.5 in a duct with e = 0.1. The amplitude B was taken as 3.5, 
and b was chosen to be 0.2. The shock path for these values is shown 
in figure (4.1). The shock starts its motion at zero velocity, acceler- 
ates, reaching a maximum velocity at t =s l/b, then gradually decelerates 
and approaches its terminal position. Over -6he time interval shown, the 
shock has traversed 91 percent of its total excursion. 

_ For an initial upstream Mach number of 1.5, the values of U and 
C are -1.36458 and 0.90972, respectively (a =» 5). The corresponding 
conditions downstream of the shock (using the approximate shock transi- 
tion relations) are U = -0.73655 and C = 1.03533. The variation of 
the downstream flow conditions along the shock path is shown in figure 
( 4 . 2 ). As the shock accelerates, the particle velocity decreases in 
magnitude rather rapidly and then increases as it gradually approaches a 
final value greater than its initial value. The sonic speed rapidly in- 
creases and then decays to a final value somewhat greater than the initial 
value. The values of the downstream flow conditions shown in figure (4.2) 
were obtained by setting a « 3 in equations (4.21). The values of the 
coefficients of equations (4.21) (cf. eqs. (4.22)) for this case are given 
in table ( 4 .I). 
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Equations (4.21) were prograjnmed for con^jutatlon on an IBM type-650 
digital computer. Tlie equations were programmed to compute the values 
of the flow variables at constant values of x, that is, at constant dis- 
tances downstream of the initial shock location. The variation of parti- 
cle velocity and sonic speed at x = -0.5, -1.0, -1.5 for the shock path 
under consideration is shown as the curves in figures (4.3), (4.4), and 
( 4 . 5 ), respectively, and is listed in table (4.2). At a1 1 cross sections, 
the history of the flow variables is similar. The particle velocity de- 
creases rapidly, reaches a minimum, and then gradually approaches a final 
value lower than its Initial value. The sonic speed undergoes a rapid 
increase, reaches a maximum, and then decreases as it slowly approaches 
a final value greater than the initial value. From this, it is apparent 
that the downstream disturbance that would cause the assumed shock motion 
consists of a compression followed by an expansion. The amplification of 
the disturbance as it moves upstream is evident in the figures. At 
X = -1.5, the magnitude of the maximum change in u is 0.0223 units or 
5.87 percent of the initial value of the variable. At x = -0.5, the 
maximum change in u is 0.0360, which is 4.55 percent of its initial 
value. Also of interest is the fact that the greater part of the change 
in the flow variables occurs in the first 5 to 6 units of time. During 
this interval, the shock wave has completed only about one- third of its 
total movement. — . 

One of the purposes of the computation of the flow field behind a 
moving shock was to determine the u,e-relatlon on lines of constant x. 
Figures (4.3) to (4.5) have been cross-plotted, and the resulting u,c- 
relations are shown in figure (4.6). At all values of x, the u,c- 
relatlon consists of an almost linear variation during the compreselon 
process. At the beginning of the expansion process, the curve folds back 
on itself and follows an almost straight path slightly displaced from that 
during the compression process. At the larger values of |xj, the portion 
of the curve representing the expansion process moves closer. to the sec- 
tion representing the compression (i.e., the "hairpin" becomes narrower). 

In section III it was assumed that, in the flow field behind the 
shock, the value of Q at a fixed x remained constant throughout the 
transient. This requires, first, that the - u, c-relatlon at constant x 
should be a straight line and, second, that the slope of the straight 
line be equal to o. In figure (4.6) the u, c-relatlon is not straight 
throughout the entire transient"; This is principally the result of the 
fact that the shock path chosen requires an e:kpansion as well, as a com- 
pression in the downstream disturbance. (Sincceeding examples will, show 
that, for shock paths which can be produced by compressive disturbances, 
the relation can be approximated by a straight line.) The portion of the 
curves of figure (4.6) representing the compression process are approxi- 
mately straight lines. For these portions of the curves, the slopes have 
been measured suad are listed below. 


5239 



5239 


NACA m 1439 


71 


X 

du/dc 

-0.5 

5.00 

-1.0 

5.12 

-1.5 

5.26 


The value of o for the. example is 5. From the table, it is seen that 
close to the initial shock location the perturbation- series solution gives 
a value of 5 for du/dc. At greater distances, the slope increases. 
Therefore, only close to the initistl location of the shock does the u,c- 
relation during the compression approximate Q = constant. Thus, for 
the region of applicability of the first-order interaction solution of 
section III, the assumption of constant Q appears to be reasonable. 

Since, for the shock to move in the prescribed fashion, the disturb- 
ance causing the motion consists of a compression followed by an expan- 
sion, the shock comes to rest sooner than it would if the disturbance 
were purely compressive in nature. If the path were such that a longer 
time interval were required for the same over-all shock displacement, the 
disturbance should approach a pure compression. To d^onstrate this, 
another case was computed for the same initial conditions and shock ex- 
cursion, but with b = 0.1, one- half the previous value. The shock path 
for this value of b is shown in figure (4.7). The time interval for 
the shock to complete 91 percent; of its travel is approximately twice 
that for b = 0.2. The history of the flow variables on the downstream 
side of the shock during its motion is shown in figure (4.8). (The coef- 
ficients of eq,s. ( 4 . 21 ) are again given by table (4.1)). Coiiparison of 
this figure with figure (4.2) (b = 0.2) shows that, although the curves 
are quite similar in shape, the magnitude of the maximum excursion in 
the flow variables is much less in the present case. This reflects the 
smaller shock velocities during the transient for the lower value of b. 

The variations of u and c with time at constant x for this 
example are shown in figures (4.9) and ( 4 .IO) for x = -0.5 and x = -2.0, 
respectively, and are listed in table (4.3). At x = -0.5, the flow vari- 
ables exhibit a sharp rise followed by a gradual approach to their termi- 
nal values. At this section the disturbance is purely compressive. Over 
80 percent of the change in the flow variables occurs in the first ten to 
twelve units of time. During this time period, the shock has moved but 
one- third of its total travel. At x = -2.0, a slight expansion follows 
the conpression process. The magnitude of this expansion is, however^ a 
much smaller proportion of the total change in the variables in this case 
than for that with b = 0.2. 
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The u, c-relation at the two values of x is shown in figure (4.11). 
At X *» -0.5, the relation is essentially a straight line over most of the 
range of the variables. The fold, characteristic of the expansion in the 
first example, does not exist. There is, however, an increase of curva- 
ture towards the very end of the transient.- At this value of x, the 
slope of the straight portion of the curve is 5.00, as it was at the same 
station for the larger value of b. At x = -2.0 (fig. (4.11b)), -where 
there is an eacpansion process of small magnitude following the compres- 
sion, the u, c-relation folds back on itself- as in the first example. 

The major portion of this curve is essentially a straight line with a 
slope of 5.48; this again indicates the increase of the- slope of the 
straight portion of the u, c-relation with distance from the initial 
shock location . 

2. Comparison with the me thod-of- characteristics solution . - To check 
the accuracy of the representation of the flow field behind a moving shock 
wave afforded by the perturbation- series solution, the flow field was com- 
puted by the method of characteristics . For this computation the condi- 
tions on the downstream side of the moving shock wave were determined from 
the approximate shock transition relations (eqs. (3.1) and (3.2)) so that 
the comparison could be made for the same initial conditions. The char- 
acteristic net for the "ejqponentlal" shock path with B = 3.5, b = 0.2, 
and an initial shock Mach number of 1.5, is shown in figure (4.12). The 
values of the flow variables at constant x were obtained by linear in- 
terpolation between the net points and are shown as the data symbols on 
figures (4.3) to (4.5). At all values of x, the variation of u an c 
with time as determined by the method of characteristics is quite similar 
to that given by the perturbation-series solution. There is an initial 
rapid compression followed by an e: 5 »ansion process. At x = -0.5, the 
two solutions are in excellent agreement during the first eight units of 
time. After this time period, the values of c given by the method of 
characteristics are greater and the values of u smaller than those of 
the perturbation-series solution. At greater distances from the InitisO. 
shock location these errors increase. The magnitude of-the differences is 
always greatest towards the end of the transient. Further, the errors 
are proportionally greater for the sonic speed than they are for the 
particle velocity. 

These errors result primarily from the representation of the under- 
lying steady flow afforded by the perturbation series rather than from 
the representation of the transient process.- It will be recalled from 
section II that the convergence and accuracy of the series solution for 
steady subsonic flow were poorest when the boundary conditions were 
specified at x = L, that is, at a high-velocity section. This is seen 
in figure (2.13), where even the third-order solution is accurate only 
for relatively short distances downstream from the boundary. This figure 
also shows that the accuracy of the representation of the sonic speed is 
much poorer than that of the particle velocity. The rather rapid decrease 



liACA M 1439 


73 


In accuracy vltli distance from an upstream boundary causes not only tbe 
errors at the beginning of the transient but also accounts for the in- 
crease of the error vith time at a given value of x. 

The coordinate axes for the inverse- shock problem are located so 
that the origin coincides "with the initial shock location. Therefore^ 
during the course of the transient the shock vave, or carrier of the 
initial data, moves upstream from the origin, and a constant x corre- 
spends to greater distances from the initial data carrier as time in- 
creases. The error in the magnitudes of the flov variables should, 
therefore increase vith time at a constant x, as observed in figures 
(4.3) to (4.5). 

To indicate the order of magnitude of the error that can be eaqo^cted 
in the original steady flov, the second-order steady-flow solution was 
conrputed for a duct with € = 0.1, with the initial data given at the 
cross section with a Mach number of 0.7. This Mach number is close to 
that behind a shock at Mach 1.5 for cr ® 5. The results are shown in 
figure (4.13) as the solid curves. The exact values axe shown as the 
dashed curves. The rapid increase of the error in c with Ixj is quite 
apparent. Prom these curves it is easily seen that a large proportion of 
the error in the flow variables observed in figures (4.3) to (4.5) is 
attributable to the Inaccuracy of the representation of the underlying 
steady flow. 

To illustrate the effect of initial shock Mach number on the results, 
an example was calculated for an "exponential" shock path with B * 3, 
b = 0.1, e s= 0.1 and an initial shock Mach number of 1.4. The shock path 

for this case is shown in figure (4.14) and the flow conditions on the 

downstream side of the shock in figure (4.15). The coefficients of equa- 
tions (4. 21) for this case are given in table (4.4). The variation of 
the flow variables at x = -0.5 and -1.0 are shown in figures (4.16) and 
(4.17), respectively, and are listed in table (4.5). The histories of 
the flow variables for .this case are very similar in nature to those given 
previously for a shock initially at Mach 1.5 moving on an exponential path 
with b = 0.1. At X = -0.5, the disturbance is entirely con^ressive with 
the major portion of the change in the flow variables occurring in the 

first ten units of time. At x = -1.0, the compression process is fol- 

lowed by a very slight escpanslon process. The amplification of the dis- 
turbance as it moves upstream is quite evident. The u, c- relations at 
the two values of x are shown in figure (4.18). Again, these curves 
are almost straight lines with that at x = -1.0 showing the characteris- 
tic fold that Indicates the expansion process. 

This case was also solved by the method of characteristics. The 
results are shown as the data synbols in figures (4.16) and (4.17). As 
in the earlier example, the representation of the particle velocity is 
more accurate than that of the sonic speed. Again, the history of the' 
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flov variables given by the method-of- characteristics solution is very 
similar to that from the series solution. The errors increase vith time 
at a fixed value of x and are greater at the larger value of x. In 
this case the errors are relatively larger -than those observed for an 
initial Mach number of 1.5. The Increase in relative error is caused by 
the fact that, for an initial shock Mach number of 1.4, the downstream 
flow is specified at a cross section with higher subsonic velocity than 
for a shock at Mach 1.5. This further diminishes the accuracy of the 
series representation of the underlying steady flow. 

3. The "cosine” path . - The examples presented thus far have all been 
for the exponential" shock path given by equation (4.23). For this path 
it was seen that the u,c-relation at constant x could be approximated 
by strai^t lines for the compression and expansion phases of the disturb- 
ance. To demonstrate that such an approximation holds for other shock 
paths, equations (4.2l) were evaluated for a path given by 

B(i - cos cct), 

f(t) = > 0, cot < Oj (4.24) 

[b, cot > «. 

This path satisfies equations (4.3) and (4.4), and the magnitude of the 
excursion of the shock is given by the coefficient B. For convenience, 
this path will be referred to as the "cosine" path. 

The case investigated was for B = 3.5 and co = </20 with an ini- 
tial shock Nfeich number of 1.5. The total shock motion and Initial shock 
Mach number are the same as those of the first- exanple for the "exponen- 
tial" path. As the coefficients of equations (4.21) are functions of the 
initial conditions only, the values given in table (4.1) apply to the 
present example. The history of- the flow variables ab x = -0.5 and 
-1.0 are shown in figures (4.19) and (4.20), respectively, and are listed 
in table (4.6). For the shock to move on the prescribed path, the dis- 
turbance consists of a compression followed by an expansion process. The 
trigonometric nature of the shock path function is refLscted in the ap- 
pearance of the variations of u and c. The variation of u and c 
have the appearance of slightly distorted, displaced cosine functions. 

The u, c-relations at constant x are shown in figure (4.21). For 
the "cosine" path, the relation again consists of an almost straight 
section followed by an abrupt fold to a second straight section, charac- 
teristic of the expansion process. Again, the slope of the straight 
portion of the u,c- curve increases with distance from the initial shock 
location. The character of the relation between the downstream flow 
variables is thus seen to be relatively independent of the nature of the 
shock path. 
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C. Discussion 

From the examples presented, it can be seen that tbe perturbation- 
series solution gives a fairly accurate picture of the history of the 
flow variables downstream of a moving shock wave. The agreement of the 
series solution with the method-of- characteristics solution is not as 
good as that observed for the unsteady-flow examples given in section II. 
However, the discrepancies were shown to arise prim ar ily from the series 
representation of the underlying steady. flow. For subsonic flow specified 
on an upstream boundary at which the flow velocity is high, the series 
representation of the steady flow is much less accurate than that for a 
flow specified at a low-speed downstream boundary. 

For an "exponential" shock path (which approximates that obtained 
in section III for the direct interaction problem) the greater part of 
the change in the flow variables was seen to occur during the early part 
of the transient. During this early part of the transient, the shock 
moved but a small part of its total travel. This is similar to the result 
obtained in section III, where it was found that, for the disturbance 
function employed, most of the shock motion occurred after the downstream 
conditions had ceased to vary. For the "exponential" path with b = 0.1, 
the variation of the flow variables at constant x (cf. fig. (4.9)) can 
be reasonably approximated by the disturbance function given by equation 
( 3 . 33 ). The largest discrepancies for such an approximation would occur 
in the region of the "knee" of the curves of figure (4.9). 

The u, c-relatlons at constant x were seen to be approximately 
straight-line relations for the "exponential" paths that did not require 
expansion processes to bring the shock to rest. For values of x close 
to the initial shock location, the u,c-relation can be approximated by 
the Q = constant relation employed in section III. At greater distances 
from the initial shock location, the slope of straight portion of the 
u,c-relatlon Increases; this iijdlcates that the constant Q assuinption 
employed in section III is reasonably accurate only for short distances 
downstream of the shock. 

For shock paths requiring an expansion process to arrest the motion 
of the shock, the u, c-relatlon at constant x is almost a straight 
line during the compression process. At the peak of the compression, 
the u,c-relation abruptly folds back on Itself and follows a slightly 
different, but again almost straight, path during the expansion process. 
This was true not only for the "exponentisQ." path but also for the "cosine" 
path; this indicates that the relationship between the dependent variables 
is not greatly influenced by the nature of the shock path. 



76 


NACA HA 1439 


To con^jute the series solution at a given value of x required 
about 7 minutes of digital computer time. Therefore, to compute the solu- 
tion at three values of x required about 21 minutes _of machine time. 

The method- of- characteristics solution required 145 minutes of machine 
time to compute the net points covering the, same range of — x and t. 

To this time must be added that required to perform the necessary inter- 
polations,. The series solution offers, therefore, a considerable saving 
in computational effort for the solution of this problem. 


Levis Flight Propulsion laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, Aug. 22, 1958 
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APPENDIX - SiaffiOIS 
A cross-sectional area of cut 

B an5)lltude 

coefficients, eq. (3.32) 

b dimensionless constant 

C Initial value of sonic speed 

C sonic speed Immediately upstream of stationary shock 

C+ characteristic of family vlth slc^e u + c 

C- characteristic of family with slope u - c 

(k) 

C\ ’ coefficient of J-th member of k-th order terms In series for sonic 
^ speed 

c sonic speed 

c sonic speed upstream of shock 

Dj^ coefficients, eq. (3.36) 

coefficients, eq. (4.14) 

K eqs. (2.46) and (2.47) 

L duct length 

M Mach number 

B eqs. (2.46) and (2.47) 

p pressure 

Q Rlemann Invariant, u - crc 

coefficients, eq. (3.26) 

coefficients, eq. (3.29) 

T time 



78 


NACA TM 1439 


t]i(k) 

O 

t 

U 


u 


U 

U 

V 

V 

M 

^3 


X 

y 

a 

P 

r 


5 

€ 

t, 

r\ 

A 

I 


coefficient of j-th. member of k-th order terms in series for time 
time 

initial value of particle velocity- 

particle velocity Immediately upstream of stationary shock 

coefficient of j-th member of k-th order terms in series for 
particle velocity 

particle velocity 

particle velocity upstream of shock 
dimensionless shock velocity 
shock -velocity 

coefficient of J-th member of k-th order terms in series for space 
coordinate 

space coordinate 

space coordinate, two-dimensional steady flo-w 

characterlstic parameter of family -with slope u - c 

characteristic parameter of family -with slope u + c 

ratio of specific heat at constant pressure to specific heat at 
constant volume 

peflod 

"smallness" parameter 
eq,. (2.23) 
eq. (2.31) 
index 


a 


parameter 

2/(r “ 1) 
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T parameter 

01 circular frequency- 

subscripts : 


s at shock location 

^ critical condition 
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TABLE 2.1. - COEEPICIEMTS OP EQUATIONS (2.58) FOR 
U = -0.37879, C = 1.08226, AND a = 5 (r = 7/s) 


\k 

j\ 

0 

1 

2 

1 

+0.47484 

-0 .03989 

-0.01455 

2 


+1.36687 

+0.47507 

3 


+0.31687 

-0.10007 

4 


+1.12250 

+0.16204 

5 



+0.46192 

6 



+0.88876 

7 



+0.02252 

8 



-0.17960 

9 



-2.19853 

10 



-0.74323 

11 



+1.82431 
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03 

to 

CM 

lO 


TABIE 3.1. - COEFFICIENTS OF EQUATIOHS (2.69) FOR 
U = -0.73179, C = 1.04541, a = 5 




■ (k) 
'’j 


\k 

1 

2 

j\ 



1 

-0.38251 

-0.09017 

2 

+5.00000 

+1.91049 

3 


-0.33714 



\k 

0 

1 

2 

1 

+0.15000 

+0.10044 

+0.02093 

2 

+0.85000 

-2.43924 

-0.70069 

3 


+0.28697 

+4.66657 

4 



+0 .08572 

5 



-0.82350 

6 



+2.88229 

7 



+0.16470 
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TABLE 3.2. - HUMERICAL INTEGRATION OP EQUATION FOR SHOCK PATH FOR 
U = -0.73179, C = 1.04541, U = -1.36458, C = 0.90972, AND X = 0.00700 


p 

First oarder 

Second order 


a 

X 

t 

a 

.X 

t 

0 

0.02629 

0.00700 

0.02236 

0.02629 

0.00700 

0.02236 

2 

2.12527 

.03421 

2.10595 

2.12536 

.03424 

2.10602 

4 

4.26928 

.07525 

4.22668 

4.26996 

.07542 

4.22724 

6 

6.42758 

.11844 

6.36060 

6.42942 

.11895 

6.36214 

8 

8.56341 

.15490 

8.47593 

8.56651 

.15576 

8.47852 

10 

10.67989 

.18570 

10.57512 

10.68403 

.18665 

10.57859 

12 

12.77972 

.21177 

12.66034 

12.78450 

.21308 

12.66436 

14 

14.86522 

.23386 

14.73349 

14.87020 

.23521 

14-.-73769 

16 

16.93839 

.25258 

16.79622 

16.94315 

.25386 

16.80022 

18 

19.00099 

.26846 

18.84997 

19.00513 

.26957 

18.85344 

20 

21.05451 

.28195 

20.89598 

21.05770 

.28280 

20.89864 

22 

23.10025 

.29340 

22.93534 

23.10223 

.29393 

22.93697 

24 

25.13932 

.30314 

24.96900 

25.13990 

.30330 

24.96943 

26 

27 .17268 

.31141 

26.99777 

27.17174 

.31118 

26.99689 

28 

29.20115 

.31845 

29.02234 

29.19863 

.31781 

29.02010 

30 

31.22545 

.32443 

31.04332 

31.22132 

.32339 

31.03969 

32 

33.24618 

.32952 

33.06123 

33.24045 

.32808 

33.05622 

34 

35.26387 

.33386 

35.07652 

35.25657 

.33202 

35.07015 

36 

37.27895 

.33754 

37.08955 

37.27015 

.33533 

37.08189 

38 

39.29181 

.34068 

39.10067 

39.28158 

.33812 

39.09178 

40 

41.30277 

.34335 

41.11016 

41.29121 

.34046 

41.10011 

42 

43.31211 

.34563 

43.11824 

43.29930 

! .34242 

43.10712 

44 

45.32008 

.34756 

45.12514 

45.30612 

.34408 

45.11301 

46 

47.32686 

.34921 

47.13101 

47.31185 

.34547 

47.11797 

48 

49 .33265 

• .35062 

49.13602 

49.31667 

.34663 

49.12215 

50 

51.33758 

.35181 

51.14029 

51.32072 

.34762 

51.12566 

52 

53.34178 

.35283 

53.14392 

53.32413 

.34844 

53.12861 

54 

55.34536 

.35370 

55.14702 

55.32699 

.34913 

55.13109 

56 

57.34840 

.35443 

57.14867 

57.32940 

.34972 

57.13317 

58 

59.35100 

.35506 

59.15191 

59.33142 

.35021 

59.13492 

60 

61.35321 

.35560 

61.15383 

61.33313 

.35062 

61.13640 

62 

63.35510 

.35605 

63.15546 

63.35456 

.35096 

63.13764 
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OIABEE 3.3. - COEEFICIEKTS OF EQUAO?IONS (2.69) FOR 
U = -0.63481, C = 1.05801, a = 5 


* 


V 







\ k 

0 

1 

2 

1 

+0.33856 

-0.11445 

-0.04173 

2 


+3.84002 

+1.69622 

3 



-11.79578 

4 



-0.09344 

5 



+1.81471 

6 



+1.20981 

7 

1 


-0.16319 


\ k 

0 

1 

2 





1 

+0.20000 

+0.08507 

+0.01716 

2 

+0.80000 

-2.26841 

-0.60600 

3 


+0.37806 

+4.28808 

4 



+0.11520 

5 



-1.07201 

6 



+2.85871 

7 



+0.21440 
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TABLE 3.4. - NUMERICAL INTEGRATION OP EQUATION FOR SHOCK PAIH FOR 
U = -0.63481, C = 1.05801, U = -1.36458, C = 0.90972, AND X = 0.83000 



First order 

Second order 


a 

X 

t 

a 

X 

t 

0 

2.69884 

0.83000 

2.22018 

2.72931 

0.82999 

2.24933 

2 

4.79665 

.86977 

4.29379 

4.77498 

.85655 

4.27834 

4 

6.93161 

.92103 

6.39908 

6.86908 

.89706 

6.34890 

6 

9.11598 

.97238 

8.55485 

9.01770 

.93749 

8.47476 

8 

11.28114 

1.01771 

10.69487 

11.14585 

.97184 

10.58363 

10 

13.42907 

1.05779 

12.82069 

13.25611 

1.00099 

12.67757 

12 

15.56156 

1.09326 

14.93369 

15.35078 

1.02574 

14.75841 

14 

17.68018 

1.12468 

17.03511 

17.43191 

1.04672 

16.82784 

16 

19.78637 

1.15253 

19.12610 

19.50134 

1.06452 

18.88736 

18 

21.88138 

1.17723 

21.20769 

21.56065 

1.07961 

20.93828 

20 

23.96637 

1.19915 

23.28079 

23.61126 

1.09240 

22.98179 

22 

26.04237 

1.21861 

25.34627 

25.65440 

1.10324 

25.01891 

24 

28.11031 

1.23589 

27.40488 

27.69113 

1.11243 

27 .05056 

26 

30.17102 

1.25125 

29.45733 

29.72238 

1.12021 

29.07750 

28 

32.22526 

1.26490 

31.50424 

31.74895 

1.12680 

31.10042 

30 

34.27369 

1.27703 

33.54617 

33.77152 

1.13238 

33.11991 

32 

36.31694 

1.28781 

35.58364 

35.79069 

1.13710 

35.13646 

34 

38.35553 

1.29740 

37.61711 

37.80696 

1.14111 

37 .15052 

36 

40.38998 

1.30593 

39.64700 

39 . 82077 

1.14450 

39 .16245 

38 

42.42071 

1.31352 

41.67368 

41.83247 

1.14737 

41.17257 

40 

44.44812 

1.32027 

43.69750 

43.84240 

1.14980 

43.18116 

42 

46.47257 

1.32627 

45.71875 

45.85081 

! 1.15186 

45.18844 

44 

48.49436 

1.33162 

47.73771 

47.85795 

1.15360 

47.19461 

46 

50.51380 

1.33637 

49.75462 

49.86399 

1.15507 

49.19984 

48 

52.53112 

1.34060 

51.76969 

51.86911 

1.15632 

51.20427 

50 

54.54656 

1.34437 

53.78313 

53.87345 

1.15738 

53.20802 

52 

56.56032 

1.34772 

55.79512 

55.87712 

1.15828 

55.21120 

54 

58.57258 

1.35070 

57.80580 

57.88023 

1.15903 

57.21389 

56 

60.58350 

1.35335 

59.81532 

59.88287 

1.15967 

59.21618 

58 

62.59324 

1.35571 

61.82380 

61.88510 

1.16022 

61.21811 

60 

64.60191 

1.35782 

63.83136 

63.88699 

1.16068 

63.21974 

62 

66.60963 

1.35969 

65.83809 

65.88859 

1.16107 

65.22113 
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lABIil 4.2. - SERIES SOLUTION FCiR PLOW FIELD Tnmrw n 
"EXPONENTIAL" SHOCK PATH 
Initial Shock Mach Number » 1.5 
B = S.5, h = 0.2, e = 0.1, a = 5 


p 


X = -0. 

5 



X = 

-1.0 



a 

t 

-u 

c 

a 

t 

-u 

c 

0 

-1.77450 

-1.48721 

0.67168 

1.04407 

-3.31655 

-2.73255 

0.61985 

1.05022 

1 

-.73179 

-.44361 

,66077 

1.04624 

-2.23927 

-1.65384 

.61071 

1.05200 

z 

.26073 

,55106 

.65302 

1.04780 

-1.22390 

-.63541 

.60419 

1.05329 

3 

1.21853 

1.51527 

.64778 

1.04887 

-.24771 

.34441 

.59975 

1.05418 

4 

2.15378 

2.46068 

.64443 

1.04958 

.70227 

1.29955 

.59692 

1.05476 

5 

3.07584 

3,39503 

,64245 

1.05002 

1.63276 

2.23926 

.59529 

1,05511 

6 

3.99107 

4.32366 

.64144 

1.05026 

2.55195 

3.17023 

.59449 

1.05529 

7 

4.90392 

5.25029 

.64109 

1.05037 

3.46S74 

4.09718 

.59420 

1.05537 

0 

5.81751 

6.17755 

.64118 

1.05039 

4.37825 

5.02342 

.59444 

1.05536 

9 

6.73401 

7.10728 

.64155 

1.05036 

5.29235 

5.95130 

.59484 

1.05531 

10 

7.65406 

8.04069 

.64208 

1.05028 

6.21002 

6 . 88238 

.59538 

1.05523 

11 

0.58100 

8.97857 

.64269 

1.05019 

7.13254 

7.81769 

.59599 

1.05513 

12 

9.51292 

9.92135 

.64334 

: 1.05009 

0.06064 

8.75781 

.59662 

1.05502 

13 

10.45084 

10.06920 

.64398 

1.04990 

0.99477 

9.70307 

.59725 

1.05492 

14 

11.39474 

11.82209 

1 .64459 

1.04988 

9.93490 

10.65^48 

.59784 

1.05482 

15 

12.34443 

12.77988 

^ .64517 

1.04970 

10.88118 

11.60894 

.59839 

1.05472 

16 

13.29964 

13.74232 

' .64670 

1.04970 

11.03313 

13.56923 

.59890 

1.05463 

17 

14,25998 

14.70908 

.64618 

1.04961 

12.79048 

13.53404 

.59936 

1.05455 

18 

15.22507 

15.67984 

.64661 

1,04954 

13.75285 

14.50304 

.59978 

1.05448 

19 

16.19447 

16.65423 

.64700 

1.04947 

14.71981 

15.47585 

.60014 

1.05442 

20 

17.16778 

17.63189 

.64734 

1.04941 

15.69092 

16.45212 

.60047 

1.05436 





p 


X = -1 

.5 



a 

t 

-u 

c 

0 

-4.72877 

-3.83893 

0.57608 

1.05470 

1 

-3,62947 

-2 .73757 

.56812 

1.05619 

2 

-2 .59981 

-1.70367 

.56243 

1.05727 

3 

-1.61501 

-.71374 

.55857 

1.05801 

4 

-.65878 

.24769 

.55610 

1.05050 

5 

.27966 

1.19117 

.55467 

1.05879 

6 

1.20404 

2.12365 

.55401 

1.05893 

7 

2,11984 

3.05046 

.55387 

1.05898 

8 

3.03222 

3.97549 

.55408 

1.05895 

9 

3.94477 

4.90146 

.55451 

1.05889 

10 

4.85995 

5.83024 

.55507 

1.05880 

11 

5.77943 

6.76306 

.55568 

1.05869 

12 

6.70423 

7.70066 

.55632 

1.05858 

13 

i 7,63492 

8.64340 

.55693 

1 1.05847 

14 

1 8.57172 

9.59142 

' .55752 

1.05836 

15 

i 9.51464 

10.54464 

.55807 

j 1.05826 

16 

: 10.46348 

11.50285 

.55850 

1 1.05017 

17 

11.41795 

12.46577 

.55903 

i 1.05809 

18 

12.37766 

13.43305 

.55944 

1,05801 

19 

13.34221 

14.40433 

.55981 

1 1.05794 

20 

14.31115 

15.37922 

.56013 

i 1.05788 


5239 



CE-12 


NACA mt 1439 


89 


TABIE 4.3. - SEEIES SOLUTIC0Sf FCR FLOW FIELD HEHIKD 
"EXPOKENTIAL" SHOCK PAIH 
Initial Shock Mach Nuniber = 1.5 


B = 3.5, h = 0.1, e = 0.1, a = 5 


p 

X = -0.5 


X = - 

-2.0 



a 

t 

-u 

c 

a 

t 

-u 

c 

0 

-1.77449 

-1.48721 

0.67168 

1.04407 

-6 .06210 

-4.85682 

0.53801 

1.05789 

2 

.23813 

.52627 

.66604 

1.04520 

-4.01271 

-2.80519 

.53436 

1.05854 

4 

2.20876 

2.50203 

.66169 

1.04609 

-2.01305 

-.80075 

.53157 

1.05905 

6 

4.15117 

4.45427 

.65841 

1.04678 

-.04556 

1.17249 

.52949 

1.05943 

8 

6.07667 

6.39225 

.65596 

1.04731 

1.89872 

3.12480 

.52800 

1.05971 

10 

7.99307 

8.32244 

.65414 

1.04771 

3.82618 

5.06339 

-52697 

1.05990 

12 

9.90582 

10.24450 

.65281 

1.04802 

5.74448 

6.99577 

-52630 

1.06003 

14 

11.81868 

12.17613 

.65183 

1.04826 

7.65907 

8.92565 

.52587 

1.06010 

16 

13.73423 

14.10994 

.65112 

1.04844 

9.57373 

10.85190 

.52562 

1.06014 

18 

15.65416 

16.03016 

.65060 

1.04858 

11.49103 

12 .78483 

.52550 

1.06016 

20 

17.57950 

17.97560 

.65023 

1.04869 

13.41267 

14.72463 

.52546 

1.06016 

22 

19.51081 

19.91099 

.64996 

1.04877 

15.33968 

16.66934 

.52548 

1.06015 

24 

21.44832 

21.86816 

.64977 

1.04883 

17 .27261 

18.60309 

.525541 

1.06013 

26 

23.39199 

23.81679; 

.64963 

1.04888 


20.55725 

. 52562 i 

1.06011 

28 

25.34159 

25.77517 

.64954 

1.04892 


22.51121 

.52571 

1.06009 

30 

27 .29682 

27.74060! 

.64948 

1.04895 


24.47308 

.52580! 

1.06006 

32 

29.25727 

29.70976 

.64944 

1.04897 

25.06418 

26.4400li 

.52590^ 

1.06004 

34 

31.22253 

31.67907 

.64942 

1.04898 


28.40870 

.52599 

1.06002 

36 

33.19213 

33.65483 

.64941 

1.04899 

28.99205 

30.38558 

.52607 

1.05999 

38 

35.165651 

35.63334 

.64941 

1.04900 

30.96258 

32.35703' 

.52616 

1.05998 

40 

37. 14270 j 

37.61138 

.64941 

1.04901 


34.33950 

.52623 

1.05996 
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OIABIE 4.4. - COEFFICIENTS OP EQUATIONS (4.21) FOR 
, U = -1.29987, C = 0.92848, U = -0.77119, C = 1.03421 
e = 0.1, (j = 5, Eg = -0.04598, Eg = 1.24223, = -0.77151 
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TABLE 4.5. - SERIES SOLUTION FOR PLOW FIELD BEHIND 
"EXPONElilTIAL” SHOCK PATH 


Initial Shock Mach Number = 1.4 


B = 3.0, h = 0.1, e = 0.1, 0=5 


3 

X = -0.5 


X = - 

-1.0 



a 

t 

-u 

c 

a 

t 

-u 

c 

0 

-1.92014 

-1.63773 

0.69799 

1.04450 

-3.53678 

-2.96193 

0.64161 

1.05140 

2 

.09526 

.37839 

.69320 

1.04545 

-1.50456 

-.92850 

.63770 

1.05215 

4 

2.06954 

2.35679 

.68953 

1.04619 

.48409 

1.06267 

.63473 

1.05273 

6 

4.01594 

4.31122 

.68680 

1.04675 

2 .44137 

3.02605 

.63252 

1.05316 

8 

5.94544 

6.25110 

.68478 

1.04717 

4.37785 

4.97192 

.63094 

1.05347 

10 

7.86566 

8.18294 

.68331 

1.04749 

6.30224 

6.90760 

.62981 

1.05369 

12 

9.78194 

10.11133 

.68224 

1.04772 

8.22072 

8.83829 

.62902 

1.05384 

14 

11.69800 

12 .03947 

.68147 

1.04789 

10.13758 

10.76764 

.62848 

1.05394 

16 

13.61640 

13.96956 

.68092 

1.04801 

12.05585 

12.69820 

.62812 

1.05400 

18 

15.53883 

15.90307 

.68052 

1.04810 

13.97754 

14.63166 

.62788 

1.05404 

20 

17 .46637 

17.84093 

.68024 

1.04816 

15.90394 

16.56914 

.62773 

1.05406 

22 

19.39958 

19.78365 

.68004 

1.04821 

17.83582 

18.51127 

.62764 

1.05407 

24 

21.33873 

21.73144 

.67989 

1.04824 

19.77353 

20.45836 

.627591 

1.05407 

26 

23.28381 

23.68431 

.67980 

1.04827 

21.71715 

22.41047 

.62758 

1.05406 

28 

25.23462 

25.64210 

.67973 

1.04828 

23.66653 

24.36749 

.62758 

1.05405 

30 

27 .19089 

27.60456 

.67969 

1.04829 

25.62145 

26.32921 

.62760 

1.05404 

32 

29.15223 

29.57137 

.67966 

1.04830 

27.58155 

28.29533 

.62762 

1.05403 

34 

31.11824 

31.54219 

.67965 

1.04831 

29.54640 

30.26549 

.62765 

1.05402 

36 

33 .08849 

33.51665 

.67964 

1.04831 

31.51560 

32.23935 

.62767 

1.05401 

38 

35.06256 

35.49439 

.67964 

1.04831 

33.48873 

34.21654 

.62770 

1.05400 

40 

37 .04004 

37 .47506 

.67964 

1.04831 

35.46538 

36.19672 

.62773 

1.05399 
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lABLE 4.6. - SERIES SOLUTION FOR PLOW FIELD BEHEND 
"COSINE" SaXK PATH 
Initial Shock Mach Nxjmber = 1.5 ■ 


B = 3.5, 0 ) = Jt/20, € = 0.1, cr = 5 


p 

X = -0.5 


X = ■ 

•1.0 



a 

t 

-u 

c 

a 

t 

-u 

c 

0 

-1.77449 

-1.48721 

0.67168 

1.04407 

-3.31656 

-2.73255 

0.61985 

1.05023 

1 

-.75764 

-.47004 

.66762 

1.05488 

-2.28658 

-1.70208 

.61646 

1.05089 

2 

.24558 

.53406 

.66334 

1.04574 

-1.27007 

-.68427 

.61288 

1.05159 

3 

1.23427 

1.52548 

.65900 

1.04661 

-.26692 

.32079 

.60923 

1.05231 

4 

2.20861 

2.50487 

.65472 

1.04749 

.72269 

1.31329 

.60564 

1.05303 

5 

3.16962 

3.47304 

.65063 

1.04833 

1.69787 

2.29364 

.60220 

1.05372 

6 

4.11850 

4.43099 

.64684 

1.04912 

2.65942 

3.26251 

.59902 

1.05436 

7 

5.05658 

5.37981 

.64344 

1.04983 

3.60841 

4.22076 

.59618 

1.05495 

8 

5.98532 

6.32066 

.64053 

1.05045 

4.54611 

5.16940 

.59375 

1.05545 

9 

6.90620 

7.25476 

.63817 

1.05096 

5.47392 

6.10955 

.59179 

1.05586 

10 

7 . 82083 

8.18342 

.63638 

1.05136 

6.39340 

7 .04246 

.59034 

1.05617 

11 

8.73087 

9.10800- 

.63522 

1.05163 

7.30619 

7.96946 

.58943 

1.05637 

12 

9.63804 

10.02992 

.63468 

1.05178 

8.21407 

8.89201 

.58905 

1.05647 

13 

10.54415 

10.95066 

.63475 

1.05180 

9.11886 

9.81162 

.58921 

1.05646 

14 

11.45103 

11.87178 

.63542 

1.05170 

10.02260 

10.72999 

.58988 

1.05636 

15 

12.36060 

12.79489 

.63665 

1.05149 

10.92722 

11.64877 

! .59103 

1.05616 

16 

13.27477 

13.72163 

.63840 

1.05117 

11.83482 

12.56975 

.59262 

11.05587 

17 

14.19551 

14.65369 

.64062 

1.05075 

12.74750 

13.49473 

.59461 

'1.05550 

18 

15.12478 

15.59276 

.64324 

1.05025 

13.66733 

14.42555 

.59694 

1.05506 

19 

16.06438 

16.54044 

.64622 

1.04968 

14.59633 

15.36399 

.59955 

1.05457 

20 

17.01610 

17.49829 

.64947 

1.04905 

15.53644- 

i 

16.31174 

.60239 

1.05403 
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Figure 2.5. - Sonic speed as a function of x. Sutsonic flow specified on 
X = O. "Epsilon duct": o = 0.1> L = 6^ a = 1. '■ 
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6 


(a - p) 
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Figure 2.6. - x as a function of (a - p) . Subsonic flow specified on x 
duct”: c c 0.1, L « 6, a = 1. 
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Figure 2.7. - P = 0 Characteristic. Subsonic flow specified on x = 0. 
"Epsilon duct": 6 = 0.1, L = 6, a = 1. ~ 
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(b) Sonic speed. 


Figure 2.8. - Concluded, 
flov specified on x = 


Flow variables as a function of x. Subsonic 
0. "Epsilon duct": e = O.Ijl. L = 6, c = 5. 
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(a) Particle velocity. 

Plguz^ 2.10. - Flow variables as a function of x. Subsonic flow specified on x » 0. "Epsilon 
duct": 0.4, L« 1.5, a r. 5. 









1 2 3 4 5 6 

X 

(a) Particle velocity. 

Figure 2.11. - Flow varial)les as a function of x. Subsonic flow spec- 
ified on X = 0. "Epsilon duct”: 0.2, L = 5, a = 5. - 
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(b) Sonic speed. 

Figure 2.11, - Concluded. Flow variables as a function of x. Subsonic 
flow specified on x = 0. "Epsilon duct”: e = 0,2 jl.L = 5, a = 5. 
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Figure 2.13. - Flow varialiles as a function of x. Subsonic flow specified 
on X = L. "Epsilon duct"; e = 0.1,L=6,ff=5. 
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(b) Sonic speed. 

Figure 2.13. - Concluded. Flow variables as a function of x. Subsonic 
flow specified on x o L. "Epsilon duct": e = 0.1,,, L = 6, ff = 5'. 
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Figure 2.14. - x as a function of (a - p) , Subsonic flow specified on x * L. "Epsilon 
duct": e = 0.1, L = 6, a = 5. 
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(a) Particle velocity. 


Figure 2.15. 
Ifled on 


- Flow variables as a function of x. Supersonic flow spec- 
: = L. "Epsilon duct": ® = 0.1^ L = 6, a « 5. 




(b) Sonic speed. 


Figure 2.15. - Concluded. Flow variables as a function of x. Supersonic 
flow specified on x = L. "Epsilon duct”: e = 0.1, L = 6, 0 = 5. 
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0 1 2 3 4 5 6 


X 

("b) Sonic speed. 

Figure 2.16. - Concluded. Flow variables as a function of x. Supersonic 
flow specified on x = 0. "Epsilon duct": e = 0.1, L = 6, cf = 5. 



Figure 2. 17- - X as a function of (a - p). Siq>ersonlc flov epecified on x 
e e 0.1, L = 6, 0=5. 
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Figure 2.22 - Variation of u and c on x =: 1.0 for ”siniple-vave-type'* disturbance 

on X 53 0; s = 0.1, c a 5. 
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(x - X) 


Figure 3.1. - Variation of flow variables upstream of a shock at 
Mach 1.5 in an "epsilon duct"; 6=0.1^ o=5. (Test of eq. 
(3.20)). 








Figure 3.3. - Variation of dovnstreani flow varla'bles in an Epsilon duct' 
Boundary Mach niimher, 0.6j ^ « 0.1, a «= 5. - 






I 


5239 





0 10 20 30 40 50 60 70 

t 


Figure 3.4. - Shock path, for a "simple- wave-type" dlsturhance specified at a downstream 
Mach ntmiber of 0.7. Initisil upstream Mach number, 1.5; ss 0.1, a a 5. 
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Figure 3.5. - Shock path for a "slmple-wave-type” dlsturhance specified at a downstream 
Mach numher of 0.6. Initial upstream Mach niimher, l-5ic=0.1,cr«=5. r 


























Pigure 4.5. - Valuation of flow rarlables at x «■ -1 
Mach number. 1.5j B =x 3.5, b =* 0.2, s a 0.1. 
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Figure 4.6. - Continued. u^c-Relatlon at constant x for 
"exponential" . shock path. Initial shock Mach number, 1.5; 
B = 3.5, b =s 0.2, e = 0.1. 
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1.044 1.045 1.046 1.047 1.048 1.049 

c 

(a) X = -0.5. 

Pigure 4.11. - u^c-Relation at constant x for 
"eajjonentlal" shock path. Initial shock Mach num- 
hcr^ 1.5j B = 3.5^ h s 0*1^ s ^ 0.1. 


Figure 4.11. - Concluded. u,c-Relatlon at 
constant x for "exponential" shock path. 
Initial shock Mach numher^ 1.5j B =• 3.5, 
h =» 0.1, c « 0.1. 
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Figure 4, 13. - Variation of flow variables with x for a subsonic steady- 
flow specified at Mach 0.7 sectlonj 6 = 0.1, u t= 5. 
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Figure - Shock path fbr ’’exporuential” function^ B « "b = 0.1; e « 0.1. 
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ngure 4.15. - Varlatlott of downstreaBi flow condltlonfl aloaag "esponential'' ohjck patb. Initial shock >Iach mniber, 
1.4; B . 3.0, h « 0.1, S =. 0.1. 
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Figure 4.17, - Variation of flow Tariables at x = -1. 
number, 1.4; B « 3.0, b = 0.1, e = 0.1. 



for "e3£i»oneiitial'* shock path. Initial shock Mach 






c 

X ""OvS* 

Figure 4-18. - Concluded. u,c-Relation at con- 
stant X for "exponential” shock path. Initial 
shock Mach xTuaiber, 1.4; B » 3-0, h « 0.1, a - O.i 

Jlgure 4.10. - u,c-Relatlon at constant x for ^exponential” 
ahock path. Initial shock Mach number, 1.4j B ^ 3.0, 
b 3“ 0.1, e = 0.1. 




Figure 4,19. - Variation of flow rariables at x = ^*5 for "coBjLne" elbocJE path- Initial shock Mach number 
1*5 j B ■ Q) =a “ 0 * 1 » 










